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ABSTRACT 


LASA  subarray  amplitude  anomalies  are  investigated  using  395  medium- 
sized events,  distributed  in  ten  azimuthally  divided  sectors.  Although  LASA 
magnitudes,  when  averaged  over  all  azimuths,  are  only  slightly  biased  relative 
to  NEIS  magnitudes,  the  amount  of  bias  varies  with  azimuth  and  subarray,  sug- 
gesting that  a simple  station  correction  for  m^  bias  is  not  adequate.  In 
addition,  fluctuations  among  LASA  subarrays  are  about  0.15  in  standard  devia- 
tion, even  when  these  magnitudes  are  calibrated  in  sectors.  Details  of  such 
fluctuations  are  explained  in  part  by  local  crustal  and  upper-mantle  hetero- 
geneities under  LASA.  The  amplitude  anomalies  are  linearly  related  to  the 
amount  of  travel-time  anomalies  iu  each  sector,  implying  that  both  effects  are 
due  to  crustal  focussing. 

Using  a fixed  effects  model,  the  authors  attempt  to  separate  the  cause 
of  m^  bias  into  sector  (azimuth)  effect,  subarray  effect,  and  subarray-sector 
interaction.  However,  even  this  detailed  modeling  could  not,  with  confidence, 
explain  m^  bias.  Using  the  reciprocity  principle,  even  with  a well  calibrated 
station,  a factor  of  two  in  uncertainty  results  for  predictions  of  a station's 
amplitude  for  an  event  only  50  km  away  from  a calibration  event  in  a region  as 
complex  as  LASA.  The  event  magnitude  uncertainty  would  probably  be  reduced  by 
network  averaging,  or  by  a non-statistical  detailed  crust  and  mantle  structure 
which  could  be  analyzed  by  ray-tracing  to  remove  source  and  receiver  effects. 
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INTRODUCTION 


Earthquake  P-wave  amplitudes  received  at  various  stations  show  a range  of 
fluctuations.  Even  when  transformed  into  a logarithmic  scale,  earthquake  mag- 
nitudes for  one  event  are  commonly  observed  to  have  standard  deviations  of 
more  than  one-half  log  unit.  Commonly  understood  causes  for  magnitude  fluctu- 
ations include:  (1)  radiation  pattern  of  the  seismic  source,  (2)  regional 
effect  of  upper-mantle  absorption,  (3)  local  crustal  structure  under  the  sta- 
tion, and  (4)  imprecise  distance-correction  terms  for  some  stations.  The 
effect  of  radiation  pattern  is  minimized  when  the  seismic  event  is  a single 
explosion  source.  Yet,  a study  of  the  LONGSHOT  explosion  (Lambert  et  al. , 

1969)  showed  that  P-wave  amplitude  scatter  of  this  explosion  had  a standard 
deviation  of  0.44  m^. 

Regional  variation  of  body  wave  magnitude  in  the  United  States  is  well 
established.  In  a recent  paper  Evernden  (1977)  gave  a concise  summary  of  past 
studies  concerning  station  magnitude  bias,  stating  that  magnitude  biases  can 
be  correlated  with  various  factors  such  as  travel-time  anomalies,  crustal 
thickness  under  the  station,  low-velocity  channel  in  the  upper  mantle,  high 
heat  flow,  and  upper-mantle  velocity  structure.  More  specifically,  Der  (1976) 
showed  that  P-wave  amplitudes  of  stations  in  the  Eastern  United  States  (EUS) 
are  generally  larger  than  those  from  stations  in  the  Western  United  States 
(WUS)  by  a factor  of  3.  Also,  he  indicated  that  most  of  the  P-wave  amplitude 
bias  can  be  attributed  to  regional  upper-mantle  absorption  related  to  the 
existence  of  a low-velocity  channel. 

Amplitude  fluctuations  of  3 to  1 are  fairly  common  when  observations  of 
an  event's  P waves  are  confined  to  a small  area.  LASA  subarray  amplitudes,  for 


Der,  Z.  A.,  1976.  On  the  existence,  magnitude  and  causes  of  broad  regional 
variations  in  body-wave  amplitudes.  SDAC-TR-76-8,  Teledyne  Geotech, 
Alexandria,  Virginia  22314. 

Evernden,  J.  F. , 1977.  Regional  bias  in  magnitude  versus  yield  measurement: 
Its  explanation  and  modes  of  evaluation.  Submitted  for  publication. 

Lambert,  D.  G. , D.  H.  von  Seggern,  S.  S.  Alexander,  and  G.  A.  Galat,  1969. 
The  LONGSHOT  experiment,  Vol.  1:  Basic  observations  and  measurements, 
SDL  Report  No.  234,  Teledyne  Geotech,  Alexandria,  Virginia  22314. 
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example,  have  long  been  known  to  register  fluctuations  of  4 to  1 for  most 
earthquakes.  Fluctuations  of  10  to  1 were  also  fairly  common  even  after  the 
LASA  reconfiguration  in  January  1974  to  a 50  km  diameter.  In  a small  instru- 
ment layout,  such  as  LASA,  the  effects  of  radiation  pattern,  regional  absorp- 
tion, and  distance  correction  term  are  nearly  identical  among  all  subarrays. 

Therefore,  amplitude  fluctuations  at  LASA  subarrays  are  strongly  affected  by 
the  local  crustal  and  upper  mantle  structure.  Note  that  these  fluctuations 
are  largely  ranaom  and,  therefore,  not  similar  to  the  significant  regional 
biases  that  Booth  et  al.  (1975),  Der  (1976),  and  North  (1977)  developed. 

These  authors  smoothed  over  the  random  fluctuations  and  established  the 

| 

regional  effects  with  statistical  confidence. 

Previous  studies  of  LASA  subarray  amplitudes  confirmed  that  local  het- 
erogeneities strongly  influence  the  P-wave  amplitudes.  Chiburis  and  Harten- 
berger  (1966)  showed  that  for  a typical  event,  signal  amplitudes  within  a 
subarray  vary  by  a factor  of  3 to  4,  and  that  signal  amplitudes  within  the 
entire  array  also  vary  by  a factor  of  3 to  4.  Klappenberger  (1967a)  studied 
the  behavior  of  LASA  signal  amplitude  variations  showing  that  amplitudes  of 
individual  sensors  and  subarrays  (phased  sums)  vary  with  an  approximately  log- 
normal distribution.  However,  his  study  of  spatial  correlation  of  amplitude 
anomalies  (1967b)  failed  to  show  a fixed  covariance  among  events.  Such  a 
result  is  expected  if  the  underlying  geology  is  not  laterally  uniform.  Indeed, 

Booth,  D.  C.,  P.  D.  Marshall,  and  J.  B.  Young,  1975.  Long  and  short  period 

amplitudes  from  earthquakes  in  the  range  0°  - 114°,  Geophys.  J.  R.  Astr. 

Soc.,  39,  523-538. 

North,  R.  G.,  1977.  Station  magntiude  bias  - its  determination,  causes,  and 
effects.  Lincoln  Laboratory,  Technical  Note  1977-24,  Lexington,  MA. 

Chiburis,  E.  F.,  and  R.  A.  Hartenberger , 1966.  LASA  signal  and  noise  ampli- 
tudes for  three  teleseismic  events,  SDL  Report  No.  151,  Teledyne  Geotech, 

Alexandria,  Virginia  22314. 

Klappenberger,  F.  A.,  1967a.  Distribution  of  short  period  P-wave  amplitude 
over  LASA,  SDL  Reprot  No.  187,  Teledyne  Geotech,  Alexandria,  Virginia, 

22314. 

Klappenberger,  F.  A.,  1967b.  Spatial  correlation  of  amplitude  anomalies. 

SDL  Report  No.  195,  Teledyne  Geotech,  Alexandria,  Virginia  22314. 
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Greenfield  and  Sheppard  (1969)  and  Iyer  (1971)  have  all  presented  interpreta- 
tions of  travel-time  anomaly  data  showing  that  the  crust  under  LASA  is  thinner 
toward  the  S30°E  direction  with  a step-like  structure  just  under  the  subarray 
AO.  However,  in  another  recent  study  (Aki  et  al.,  1976),  three-dimensional 
modeling  was  used.  Results  of  this  modeling  suggested  that  such  changes  in 
crustal  thickness  are  likely  to  be  overestimated  because  the  velocity  struc- 
ture of  layers  deduced  from  such  modeling  indicates  that  both  the  crust  and 
the  lower  lithosphere  appear  to  share  the  same  anomaly  pattern  as  the  proposed 
Moho  topography. 

Whether  or  not  crustal  thickness  changes  abruptly,  underlying  crustal 
and  upper-mantle  structure  both  strongly  influence  travel-time  and  amplitude 
anomalies.  The  shallow  structures  cause  local  variations  in  travel  time  and 
amplitude  observed  at  nearby  stations,  and  factors  lying  in  the  deeper  litho- 
sphere may  cause  a regional  effect. 

This  report,  using  LASA  subarray  data,  examines  whether  seismic  receiver 
amplitude  anomalies  and  magnitude  bias  can  be  explained  by  regional  and  local 
effects.  The  study  is  divided  into  three  parts.  First,  subarray  beam  ampli- 
tudes and  array  beam  amplitudes  are  compared  and  systematic  errors  due  to  the 
array's  operational  method  of  magnitude  estimation  are  evaluated.  Second, 
subarray  m^'s  are  compared  with  the  NEIS  m^  to  evaluate  azimuthal  variations 
of  Am^  in  terms  of  regional  and  local  effects.  Third,  a statistical  model  is 
used  to  describe  LASA  Am^  as  a means  of  evaluating  the  significance  of  local 
and  regional  effects  stemming  from  crustal  and  upper-mantle  structure  and  to 
evaluate  the  precision  of  m,  recorded  at  a single  site. 


Greenfield,  R.  J.,  and  R.  M.  Sheppard,  1969.  The  Moho  depth  variations  under 
LASA  and  their  effects  of  dT/dA  measurements,  Bull.  Seism.  Soc.  Am., 

59,  409-420. 

Iyer,  H.  M. , 1971.  Variation  of  apparent  velocity  of  teleseismic  P-waves 

across  the  Large-Aperture  Seismic  Array,  Montana,  J,  Geophys.  Res.,  76, 
(35),  8554-8567. 

Aki,  K. , A.  Christof ferson,  and  E.  S.  Husebye,  1976.  Three  dimensional  seismic 
structure  of  the  lithosphere  under  Montana  LASA,  Bull.  Seism.  Soc.  Am., 

66,  501-524. 
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DATA 


Events  for  this  study  were  selected  from  the  LASA  1974  event  summary  by 
collating  with  the  NEIS  bulletin.  The  year  1974  was  selected  to  ensure  rele- 
vance of  the  study  to  present  systems  because  in  late  1973  LASA  was  recon- 
figured to  a smaller  size.  Only  those  events  which  collated  with  the  NEIS 
bulletin  in  the  m^  magnitude  range  of  5.0  to  5.5  were  used.  This  m^  range 
limitation  restricted  LASA-recorded  amplitudes  of  the  selected  events  to 
between  10  and  100  millimicrons.  Since  LASA  instruments  are  fixed  to  satu- 
rate (clip)  at  350  millimicrons,  this  upper  and  lower  limit  forced  event 
select'd  from  the  medium  magnitude  range,  thus  avoiding  noise  contamination 
of  low-amplitude  signals  and  non-linear  instrument  response  of  high-ampj ltude 
signals.  During  the  analyses,  some  events  were  eliminated  because  of  data 
problems  or  because  the  signal  was  mixed  with  the  coda  of  another  event,  etc., 
yielding  a final  total  of  395  events  in  this  study. 

These  events  were  then  grouped  into  10  sectors  based  upon  back  azimuths 
of  the  events.  The  azimuthal  range  of  each  sector  was  adjusted  so  that  the 
sectors  were  easily  characterized  in  geophysical  terms.  Thus,  the  number  of 
events  in  each  sector  ranges  from  6 to  151.  Table  I lists  the  azimuthal 
range  of  each  sector  and  the  approximate  seismic  regions  each  sector  covers. 
These  data  are  also  shown  in  Figure  1. 

For  each  event  selected,  maximum  peak-to-peak  amplitudes  and  periods  of 
the  P arrival  within  3 seconds  of  the  onset  on  the  unfiltered  array  beam  and 
subarray  beams  from  the  A,  B,  C,  and  D rings  were  measured  from  waveform 
records  in  the  Seismic  Data  Analysis  Center  (SDAC)  event  library.  Chiburis 
and  Hartenberger  (1966)  showed  that  subarray  beams  resulted  in  about  2 db  of 
signal  degradation.  However,  noise  reduction  in  subarray  beams  is  approxi- 
mately 6 db  with  respect  to  individual  sensors. 

Magnitude  biases  are  computed  by  subtracting  NEIS  m^  from  calculated  LASA 
m^.  NEIS  reported  magnitudes  are  rounded,  not  truncated,  to  the  nearest  0.1 
magnitude  unit.  There  is  some  question  that  "small"  NEIS  m^  can  represent 
the  "true"  magnitude  because  of  the  small  number  of  stations  (say,  less  than 
10)  reporting  amplitudes  for  the  event.  However,  magnitudes  of  moderately 
large  events,  such  as  those  in  this  study’s  data  base,  are  generally  close  to 
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TABLE  I. 


s 

Division  of  Data  Base 

Into  10  Sectors 

- 

Sector 

Azimuthal  Range 

General  Areas 

9 

1. 

350.1 

- 

25.0 

Alma  Ata  to  Caspian  Sea 

2. 

25.1 

- 

65.0 

Turkey , Greece , Spain 

3. 

65.1 

- 

110.0 

North  and  South  Atlantic  Ridge 

4. 

110.1 

- 

130.0 

Puerto  Rico 

5. 

130.1 

- 

160.0 

South  of  Mexico  to  South  of  Chile 

6. 

160.1 

- 

210.0 

Mexico,  Easter  Islands 

7. 

210.1 

- 

250.0 

Tuamotu,  Samoa,  Fiji 

8. 

250.1 

- 

290.0 

New  Hebrides  to  Solomons 

9. 

290.1 

- 

325.0 

Aleutians,  Kurils,  Japan, 

Marianas 

t 

10. 

325.1 

— 

350.0 

Tibet,  India 
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true  magnitudes  because  many  stations  (say,  more  than  10)  report  to  NEIS  and, 
the  result  of  averaging  many  stations  is  to  closely  approximate  "true"  magni- 
tudes (Husebye  et  al.,  1974).  Since  a number  of  events  are  used  and  averaged 
in  each  sector,  the  effect  of  the  NEIS  bias  is  further  reduced.  The  fact  that 
some  NEIS  regional  bias  may  remain  is  unimportant  to  most  conclusions  of  this 
study.  Exceptions  to  this  statement  are  taken  up  as  they  arise. 

Improper  use  of  distance  correction  and  use  of  different  distance  cor- 
rection tables  may  introduce  still  more  errors.  Body  wave  magnitudes  are  com- 
puted by  the  formula  m^  = log(A*C/T)  + B(A,h),  where  A is  the  maximum  peak-to- 
peak  amplitude,  C is  the  instrument  correction  factor,  and  T is  the  period  of 
the  signal.  The  B(A,h)  factor  is  a distance  magnitude  correction  factor. 

There  are  two  types  of  B-factors  available:  Veith  and  Clawson's  P-factors 
(1972),  and  Gutenberg  and  Richter's  values  known  as  B-(or  Q-)  factors  for 
zero  to  peak  amplitudes  (converted  to  peak  to  peak  by  Cannon,  1967).  In  LASA 
operation  all  magnitudes  are  computed  with  Veith  and  Clawson's  P-factors 
(1972).  The  difference  between  P-tables  and  Gutenberg  and  Richter's  B-tables 
are  as  much  as  0.2  magnitude  unit  for  surface  focus  and  they  can  be  even 
greater  for  deeper  events.  In  addition,  two  more  computational  errors  are 
involved  in  original  LASA  m^  computations.  The  first  is  due  to  the  original 
LASA  event  location  which,  on  the  average,  is  in  error  by  about  1.5  degrees. 
Still,  this  e».ror  would  result  in  only  a small  change  in  the  distance-cor- 
rection value':.  The  second  error  is  a result  of  LASA's  failure  to  estimate 
proper  depth,  so  a wrong  correction  term  was  used  in  the  m^  computation.  To 
avoid  these  errors,  array  and  subarray  m^  have  been  recomputed  from  the  auth- 
ors's own  amplitude  (peak-to-peak)  and  period  measurements  with  the  NEIS  epi- 
center and  depth  and  with  Gutenberg-Richter  B-tables,  as  in  NEIS  magnitude 
computations. 

Husebye,  E.  S.,  A.  Dahle,  and  K.  A.  Berteussen,  1974.  Bias  analysis  of 

NORSAR-  and  ISC-  reported  seismic  event  m,  magnitudes.  J.  Geophys.  Res., 

_79,  (20),  2967-2978. 

Veith,  K.  F.  and  G.  E.  Clawson,  1972.  Magnitude  from  short-period  data,  Bull. 

Seism.  Soc.  Am.,  62,  435-452. 


Cannon,  H.  J. , 1967.  HYP01,  Technical  Report  56,  106,  Teledyne  Geotech 
Garland,  Texas. 


LASA  MAGNITUDE  FLUCTUATION  STATISTICS 


Systematic  Magnitude  Bias  Errors 

Calculating  array  magnitude  requires  first  forming  the  array  beam  (phased 
sum),  and  then  calculating  magnitude  from  it.  As  Husebye  et  al.  (1974)  dis- 
cussed, this  procedure  results  in  three  types  of  biasing  errors.  The  first 
error  results  from  signal  losses  associated  with  array  beamforming,  caused  by 
improper  travel-time  corrections,  signal  waveform  difference  among  subarrays, 
and  increased  noise  interference  on  subarrays  compared  to  the  full  beam. 
Although  beamforming  loss,  dA^ogg,  is  the  estimate  of  the  reduction  of  array 
beam  amplitude  from  the  average  subarray  amplitude,  it  does  not  comprise  the 
magnitude  bias,  because  magnitudes  are  computed  with  amplitudes  and  periods. 
The  second  measure  of  error,  dm^ogg,  is>  therefore,  the  measure  of  deviations 
of  the  array  beam  A/T  from  the  average  subarray  A/T.  The  third  type  of  non- 
random  error,  dmgkew,  is  the  difference  between  log(average  A/T)  and  average* 
log (A/T). 

Because  this  study  focuses  on  measurement  of  magnitude  bias  at  LASA,  all 
bias  computations  are  taken  with  the  following  formula: 

. LASA  NEIS 

Am . m - m /i\ 

ink  13k  (1) 

for  i-th  subarray,  j-th  sector,  k-th  event.  Using  this  fort.ula,  LASA  m^  is 
greater  than  NEIS  m^  if  Am^  is  positive,  and  negative  if  vice-versa.  Simi- 
larly, dA^ogs  and  dm^ogg  are  positive  when  array  amplitudes  (magnitudes)  are 
higher  than  subarray  averages.  Therefore,  the  sign  of  these  measurements 
are  the  reverse  of  Husebye  et  al.s'  (1974).  Finally,  the  total  non-random 
m^  bias,  dm^^,  is  the  sum  of  dm^oss  and  dmg^ew  measured  in  each  sector. 

The  following  formulas  express  the  relation  of  non-random  errors  measured 


at  LASA: 


^ias  = ^loss  + ^skew 


= log  A_  - log  - L 


where  Ag  is  the  maximum  zero-to-peak  amplitude  of  the  LASA  array  beam,  and  a^^ 
is  the  maximum  zero-to-peak  amplitude  of  the  i-th  subarray  for  that  event. 

Cg  and  c^  are  instrument  response  correction  factors,  and  Tfi  and  are  periods 
of  array  beam  (B)  and  subarray  beam  (i)  signals.  By  combining  (4)  and  (5), 


the  total  di 


"^ias 

dmbl 


in  (2)  can  be  simplified  as, 


as 


log  /_ B 


(6) 


which  could  have  been  directly  written.  Errors  with  these  formulas  are  com- 
puted for  each  independent  event  and  averaged  over  the  events  in  each  sector; 
Table  II  shows  the  results. 

Unlike  Husebye  et  al.'s  (1974)  evaluation  of  NORSAR,  where  each  region 
showed  dA^ogg  < 0,  dA^Qgg  for  LASA  shows  six  gains  and  four  losses  in  the  ten 
sectors.  This  means  that  Ag  is  slightly  higher  than  the  average  of  a_^,  the 
distribution  of  a^  is  perhaps  slightly  skewed,  and  that  the  loss  of  LASA  array 
beamforming  is  very  small.  Similarly,  dm^ogg  was  such  that  array  magnitudes 
were  biased  higher  than  the  subarray  average.  Since  dm^^  is  all  positive, 
the  total  magnitude  bias  is  all  positive,  except  in  sector  2.  Averaging  all 
sectors  the  LASA  magnitude  bias  is  0.038  m^  units,  a figure  demonstrating 
that  systematic  error  at  LASA  is  quite  small. 

Standard  Deviations  of  and  Subarray  m^ 

The  magnitude  biases  given  in  Table  II  are  computed  by  adding  two  system- 
atic losses,  dm^ogs  and  dm^^  with  formulas  (4)  and  (5).  Alternatively, 
dmg^as  and  deviations  of  dnig^ag  can  both  be  calculated  at  the  same  time  using 
formula  (6).  Thus,  simplify  formula  (6)  with  the  following  notations: 
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TABLE  II 


LASA  array  beam  losses  and  magnitude  biases 
for  each  sector  calculated  with  equations  3,  4,  5,  and  6. 


Sector 

Azimuth 

range 

No.  of 

events 

Beam  Loss 

Magnitude  Bias  1 

dA. 

loss 

dm, 

loss 

dm 

skew 

dm  . 
Dias 

1 

350.1-25.0 

35 

0.089 

-0.019 

0.022 

0.003 

2 

25.1-65.0 

22 

-0.082 

-0.048 

0.032 

-0.016 

3 

65.1-110.0 

16 

0.153 

• 

0.038 

0.034 

0.072 

4 

110.1-130.0 

6 

-0.091 

-0.006 

0.020 

0.014 

5 

130.1-160.0 

97 

-0.015 

0.020 

0.032 

0.052 

6 

160.1-210.0 

20 

0.124 

0.078 

0.023 

0.101 

7 

210.1-250.0 

24 

0.059 

0.044 

0.012 

0.056 

8 

250.1-290.0 

17 

0.069 

0.035 

0.013 

0.048 

9 

290.1-325.0 

151 

-0.072 

-0.015 

0.028 

0.013 

10 

325.1-350.0 

7 

0.075 

-0.011 

0.047 

0.036 

AVERAGE 

0.031 

0.012 

0.026 

0.038 

L 


dm, 


Bjk 


m„ ..  - m .. 
Bj  k . j k 


(6a) 


where 

dm  = magnitude  bias 
m = magnitude 
B = the  array  beam 
j = a subscript  for  the  sector 
k = a subscript  for  the  event 

. = a subscript  for  the  averaged  value  over  subarrays 
Each  independent  event  is  associated  with  a value  of  Thus,  by  aver- 

aging over  the  entire  K events  in  sector  J,  the  average  dm^^  and  error  is 
obtained,  that  is  dm„  .,  = dnu  . + e. 


lmBjk  = dmBj.  + eBjk 


(7) 


The  term  dm^  is  the  mean  magnitude  bias  for  the  j-th  sector,  and  it  is  equi- 


valent to  dm,  . shown  in  Table  II;  its  standard  deviation  a measures 
d las  Bi • 

the  width  of  the  distributions. 


Similarly,  evaluating  how  subarray  m^'s  fluctuate  from  the  mean  can  be 
done  using  the  same  subscript  notation,  let 


m.  = m ..  + e. 
ljk  ,jk  13k 


(8) 


where  is  the  magnitude  for  the  i-th  subarray,  j-th  sector,  and  k-th 


event;  and  m is  the  average  magnitude  of  all  subarrays.  For  a particular 
• 1 k 


jk-th  event,  subarray  magnitude  fluctuates  from  the  mean  with  a standard 
deviat 
yields 


deviation  of  o ,,  Averaging  over  the  entire  k events  in  the  j-th  sector 

• J K.  • 


a ..  = a . 

• jk  .3. 


+ e 

•3k 


(9) 


The  term  a . is  the  mean  standard  deviation  of  the  subarray  magnitude  fluctu- 
*J  • 


ations  in  the  j-th  sector,  and  0(0  ) represents  the  standard  deviation  of 

• j • 


the  term  e . , shown  in  (9).  Results  of  these  computations  using  (7)  and  (9) 
. j k 


are  shown  in  Table  III. 


Magnitude  biases  in  the  first  column  of  Table  III  are  almost  the  same 
as  those  in  the  last  column  of  Table  II.  Although  LASA  m^  bias  in  column  1 


are  small  (-0.031  raL  over  all  events-)  standard  deviations  are  in  the  range 
b 


1 
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TABLE  III 

Measurements  of  magnitude  bias,  and  measurements 
of  subarray  magnitude  fluctuations 


Sector 

Magnitude  Bias 

Subarray 

m Fluctuations 
b 

°Bj. 

a . 

•3. 

a(a  . ) 

• j • 

i 

0.003 

0.134 

0.144 

0.035 

2 

-0.017 

0.185 

0.177 

0.030 

3 

0.072 

0.115 

0.166 

0.057 

4 

0.014 

0.048 

0.133 

0.029 

5 

0.051 

0.092 

0.160 

0.055 

6 

0.101 

0.109 

0.134 

0.046 

7 

0.055 

0.125 

0.106 

0.016 

8 

0.047 

0.063 

0.108 

0.016 

9 

0.013 

0.151 

0.157 

0.027 

10 



0.036 

0.119 

0.213 

0.026 

of  0.1  which  means  that  even  if  LASA  biases  are  calibrated  in  each 

sector  with  dm„ . , the  accuracy  of  the  estimates  is  limited  by  the  value  of 
Bj . 

a in  column  2. 

Bj . 

Measurements  of  subarray  magnitude  fluctuations,  o . , are  in  the  range 

• J • 

of  0.11  to  0.21  m^,  with  average  0.15  m^.  Also,  the  standard  deviation 
a (a  ) is  very  small  in  each  sector,  making  subarray  m^  fluctuations  consist- 
ently in  the  range  of  0.15  m^.  The  range  of  subarray  fluctuations  — from 
smallest  to  largest  subarray  would  include  14  of  the  16  subarrays  or  - 90% 
of  the  samples.  The  extreme  samples  would  lie  at  s ±1.64o,  so  if  o = 0.15 
then  the  extreme  subarrays  for  each  event  would  be  expected  to  be  0.49  mag- 
nitude units  apart  about  a factor  of  3. 


Biases  Compared  to  NEIS  m^ 

For  events  selected  in  a sector,  the  magnitude  bias  Am^  is  computed  using 
(1)  and  then  averaged  over  the  number  of  events  in  the  sector.  The  result 
includes  magnitude  biases  for  each  subarray  and  sector,  as  well  as  associated 
standard  deviations..  These  values  are  shown  in  Table  IV.  The  last  column 


buDarray  m^ 


Magnitude  biases  may,  of  course,  be  due  in  part  to  NEIS  regional  bias. 


The  standard  deviations  of  m^  biases  shown  in  this  table  are  generally 
in  the  range  of  0.30  to  0.35.  Standard  deviations  of  all  subarrays  are  lower 
in  sectors  4 and  10.  However,  this  result  is  suspect  because  only  six  or 
seven  events,  respectively,  were  analyzed  in  each  sector.  With  the  exception 
of  these  two  sectors,  the  standard  deviations  are  stable  for  all  subarrays  and 
sectors.  Thus,  while  subarray  magnitudes  can  be  azimuthally  calibrated,  the 
best  possible  improvement  are  magnitude  estimates  with  standard  deviations  of 
0.3.  This  table  also  shows  that  the  standard  deviations  of  the  full  array 
(row  labeled  LASA)  are  in  the  same  range  as  the  subarrays,  indicating  that  any 
signal  estimate,  array  or  subarray,  can  be  azimuthally  calibrated  to  behave 
like  a single  station  while  the  standard  deviations  remain  the  same.  If  the 
NEIS  regional  bias  varied  from  event  to  event,  this  would  contribute  to  the 
variance.  However  since  5.0  <_  m^  5.5  the  bias  would  not  be  expected  to  vary 
much  because  the  same  suite  of  station  would  be  expected  to  detect  all  events. 
Earthquake  radiation  pattern  variations  are  presumably  responsible  for  much  of 
this  variance.  (Recall  that  the  a of  subarray  m^  relative  to  the  mean  LASA 
beam  m^  is  0.15.) 


Uniform  Azimuthal  Variations  of  Subarray  Magnitude  Biases  — a Regional  Effect 


Evernden  (1977)  suggested  that  P-wave  amplitudes  received  at  stations  in 


a broad  region  can  be  calibrated  with  a regional  correction.  This  observa- 
tion suggests  that  one  correction  value,  such  as  Aim  in  Table  IV,  applies  to 


events  from  all  directions.  Although  the  array  Aim  is  small  and  near  the 


value  of  0.1  in  this  table,  azimuthal  variations  are  much  greater  in  all  sub- 
arrays, suggesting  that  a simple  correction  is  not  adequate  for  any  particu- 
lar subarray.  Figure  2 shows  the  plot  of  Aim  against  azimuth  for  all  sub- 
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SECTOR 


Figure  2.  Azimuthal  variations  of  LASA  subarray  magnitude  bias.  Amb  of  all 
subarrays  follow  general  trends  indicating  such  variations  are  the 
effects  of  the  deeper  structure  which  affects  all  subarrays. 


Mean  magnitude  bias  and  standard  deviations  of 
LASA  subarrays  in  10  divided  sectors. 
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arrays,  indicating  that  azimuthal  variations  of  magnitude  bias  are  acting 
uniformly  over  most  subarrays,  although  some  variations  in  each  sector  are 
either  greater  or  less  than  the  average  magnitude  bias. 

If  factors  causing  m^  bias  in  subarrays  may  be  separated  into  regional 
and  local  structural  changes,  then  the  appearance  of  Am^  variations  acting 
uniformly  in  Figure  2 would  show  regional  variation  and  the  plot  of  Am^  with- 
in each  sector  (Figure  3)  would  reflect  local  structural  changes.  Based  upon 
this  assumption,  the  regional  cause  is  deeper  than  the  local  cause  and,  there- 
fore, the  regional  cause  may  be  the  P-wave  velocity  and  Q variations  of  the 
upper  mantle  (Evernden,  1977).  Consider  that  LASA  is  located  at  the  edge  of 
the  Western  United  States  structure  and  that  Q might  be  higher  in  the  north 
and  northeast  of  LASA  and  lower  in  the  Southwest  and  West.  Consequently, 
events  at  LASA  from  sectors  1 and  2 might  show  slightly  positive  magnitude 
bias  and  events  in  sectors  7 and  8 might  show  strongly  negative  bias  as  is 
observed.  This  is,  of  course,  only  proof  of  consistency,  and  other  causes 
such  as  source  regional  NEIS  bias,  or  source  radiation  pattern,  are  possible. 

Contoured  Patterns  of  Subarray  Sector  Magnitude  Biases  — a Local  Effect 

Figure  3 shows  contour  maps  of  average  magnitude  bias  by  sector;  the 
arrow  in  the  figure  shows  direction  of  the  incoming  signal.  Although  small, 
these  Am^  contours  show  a close  resemblance  to  similar  contours  made  of  tra- 
vel-time anomalies  (Iyer,  1971).  The  most  conspicious  pattern  in  the  contour 
maps  of  Am^  is  the  N60°E  trend,  which  is  in  good  agreement  with  the  contoured 
travel-time  anomalies  presented  by  Aki  et  al.,  (1976),  Iyer,  (1971)  and 
Greenfield  and  Sheppard  (1969). 

Interpretations  of  the  N60°E  trend  of  travel-time  anomaly  data  is  some- 
what uncertain.  Greenfield  and  Sheppard,  as  well  as  Iyer,  suggested  this 
trend  indicated  that  the  Moho  boundary  is  sloping  upward  from  under  AO  toward 
the  southeast.  However,  Aki  et  al.  (1967)  and  Capon  (1974)  suggested  that 
perhaps  such  shallowing  of  the  Moho  is  exaggerated.  While  we  can  offer  no 
conclusive  evidence,  our  interpretation  of  contour  lines  in  Figure  3 is 

Capon,  J.,  1974.  Characterization  of  crust  and  upper  mantle  structure  under 

LASA  as  a random  medium,  Bull.  Seism.  Soc.  Am.,  64,  235-266. 
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Figure  3.  Contour  maps  of  average  m^  bias  in  10  divided  sectors.  The  vari- 
ations in  contour  lines  suggest  the  influence  of  variations  in 
local  crust  upper-mantle  structure  (located  at  the  end  of  the 
report  as  a foldout). 
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consistent  with  Greenfield's  structure. 


Figures  4a,  4b,  and  4c,  show  NW-SE  cross  sections  of  LASA  cutting  from 
D4  to  D2.  The  depth  of  the  Moho  was  taken  from  the  crustal  structure  in 
Figure  4 of  Iyer's  (1971)  article.  The  18°  updip  of  the  Moho  boundary,  begin- 
ning at  AO  and  moving  toward  the  SE  is  in  agreement  with  the  model  that  Green- 
field and  Sheppard  (1969)  proposed,  except  that  in  their  paper  the  Moho  appears 
deeper  than  it  was  placed  by  Iyer.  A simple  velocity  model  of  6 km/sec  for 
the  crust  and  8 km/sec  for  the  mantle  is  used  in  this  study  to  compute  the 
incident  wave  refractions  at  the  Moho  boundary. 

In  Figure  4a,  where  seismic  rays  are  emerging  in  parallel  with  the  trend 
of  the  thinning  crust  under  LASA,  the  rays  are  slightly  convergent.  The  con- 
vergence appears  between  B2  and  AO,  which  agrees  with  the  contour  lines  in 
Figure  3 for  Sectors  2 and  7.  Note  that  this  interpretation  is  based  upon 
a simple  model  that  previous  workers  suggested.  Contour-line  features  appear 
to  be  simpler  in  the  NW  corner  (on  the  flat  crust  side)  and  more  complex 
toward  the  SE  (thinning  crust  side). 


Figure  4b  shows  the  case  corresponding  to  events  in  Sector  5,  with  the 
average  angle  of  incidence  of  seismic  waves  of  25°  (corresponds  to  events  from 
Chile-Boliva  border  area) ; waves  will  be  divergent  at  the  surface  because 
they  will  encounter  the  updip  Moho  structure.  The  divergent  feature  is  rather 
simple  in  this  case,  the  effect  of  the  varying  crustal  thickness  appears  on 
the  contour  lines  as  the  broadening  of  line  spacing.  This  feature  will  begin 
as  a N60°  trend  lying  between  D4  and  B4  and  continue  toward  the  SE.  Contour 
lines  in  Sectors  4 and  5 agree  well  with  this  explanation. 


Figure  4c  shows  waves  emerging  from  the  NW  corresponding  to  Sector  10 
contour  lines.  An  emergence  angle  of  16°  was  used  which  corresponded  to  events 
from  the  Tashkent  area.  The  upgoing  rays  are  convergent  at  the  surface  near 
B2.  The  waveforms  are  complex  because  converging  rays  and  contour  lines  are 
narrowly  spaced  with  a N60%E  trend.  Such  features  agree  well  with  Sector  10 
contour  lines  in  Figure  3. 


In  addition  to  the  NE-SW  trend  of  these  contour  lines  Figure  3 also  shows 
some  local,  small-scale  magnitude  highs  and  lows.  These  small  features  may 
be  the  shadow  of  local  heterogeneities  such  as  small  "bumps"  in  the  Moho  boun- 
dary or  small  local  low  Q zones.  As  the  azimuth  of  emerging  waves  changes 
from  one  sector  to  another,  these  features  show  corresponding  responses  on  the 
surface.  Judging  from  the  relative  movements  of  these  small  feature  against 
the  NE-SW  trend,  small  irregularities  are  likely  to  be  located  either  directly 
above  or  below  the  Moho  boundary. 

To  investigate  the  relation  between  amplitude  and  travel-time  anomalies, 
some  values  of  travel-time  anomalies,  At,  were  chosen  from  Chiburis’  and 
Ahner's  (1973)  work.  Selection  was  made  based  upon  the  location  of  events 
shown  in  Figure  1;  then  one  location  was  chosen  per  sector  where  most  events 
occurred  in  that  sector.  Although  selecting  one  representative  location  in 
a sector  was  difficult,  the  choices  may  be  adequate  for  first-order  estimates 
of  amplitude  and  travel-time  anomaly  relations.  After  determining  a repre- 
sentative location,  travel-time  anomalies  of  that  location  where  taken  from 
the  Chiburis-Ahner  travel-time  correction  table.  This  method  was  simpler 
than  the  more  precise  method  of  relating  travel-time  corrections  of  each 
subarray,  sector,  and  event,  At^^*  with  Arn^^. 

Figure  5 shows  the  plot  of  Am^  against  At  for  each  sector,  where  Am^'s 
are  values  shown  in  Table  IV,  and  At's  are  Chiburis-Ahner  correction  values. 
Although  some  wide  scattering  exists,  perhaps  a result  of  selecting  only  one 
location  per  sector,  amplitude  and  travel-time  anomalies  are  linearly  and 
clearly  related;  that  is,  positive  Am^  is  related  to  positive  At.  Note  that 
this  is  opposite  to  the  relation  discussed  in  Evernden's  (1977)  article, 
where  he  suggested  that  high  attenuation  (negative  Am^)  and  positive  (late) 

At  correlated  in  the  case  of  broad  regional  m^  anomalies.  The  linear  rela- 
tion in  Figure  5 suggests  that  the  LASA  subarray  effects,  which  are  inter- 
preted to  be  shallower  than  the  regional  effect,  are  generally  consistent 


Chiburis,  E.  F. , and  R.  0.  Ahner,  1973.  LASA  regional  travel-time  corrections 
and  associated  modes,  SDAC-TR-73-6,  Teledyne  Geotech,  Alexandria,  Virginia, 
22314. 
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Figure  5.  The  relation  of  magnitude  bias  and  travel-time  anomalies  at  LASA 
subarrays  in  each  sector. 
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with  ray  theory.  For  example,  if  a small  body  of  high  velocity,  shaped  as 
as  a lens,  exists  in  the  upper  mantle,  then  seismic  rays  through  this  lens  will 
be  divergent  and  the  ray  arrives  earlier  (negative  At)  than  normal  to  the 
receiver  above  the  lens  and  is  associated  with  smaller  (divergent)  amplitude. 
Conversely,  for  low-velocity  lens  seismic  rays  will  be  convergent  and  ray 
arrives  later  (positive  At)  with  positive  Am^. 

Because  the  low  Q zone  is  deep,  a negative  Am^  effect  will  be  uniform 
for  all  subarrays.  The  net  result  is  a shift  of  the  coordinate  origin  down- 
ward (Since  At  are  relative  to  AO  there  is  no  effect  on  At.)  In  this  context 
Figure  5 (and  Figure  2)  shows  that  Sectors  6,  7,  and  8 are  associated  with 
coordinate  shifts,  and  inspection  of  other  sectors  demonstrates  that  while 
Sector  1,  2,  and  5 have  positive  coordinate  shifts.  Sectors  3 and  4 need  no 
origin  shift.  This  overall  situation  is  consistent  with  the  theory  that  a 
low  Q zone  exists  in  the  WUS  and  that  LASA  is  situated  just  at  the  edge  of 
this  zone,  so  that  events  in  Sectors  1 and  2 seem  to  show  the  effect  of  waves 
travelling  through  Canadian-shield  type  mantle  and  events  from  the  SW,  W and 
NW  (Sectors  6,  7,  and  8)  show  the  effects  of  waves  travelling  through  the  high 
attenuation  zone  of  the  WUS  mantle. 

In  summary,  a comparison  of  LASA  and  NEIS  magnitudes  suggests  that  all 
station  magnitude  biases  should  be  calibrated  azimuthally  rather  than  with 
a single-value  bias  term.  The  data  also  shows  that  magnitude  biases  are 
affected  by  regional  upper-mantle  attenuations,  as  well  as  by  local  irregular 
structures  located  within  the  upper  100  km  of  the  earth’s  crust,  and  that 
linear  positive  relation  of  amplitude  and  travel-time  anomalies  suggests  that 
focussing  due  to  local  crustal  irregularities  is  probably  the  main  cause  of 
magnitude  anomalies  within  the  array.  Finally,  it  was  revealed  that  a stan- 
dard deviation  of  0.3  m^,  with  respect  to  NEIS,  exists  even  after  providing 
a good  local  and  azimuthal  correction;  this  is  presumably  due  to  earthquake 
radiation  patterns  and  other  source  effects  which  "focus"  a varying  propor- 
tion of  energy  toward  LASA. 
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FORMAL  ANALYSIS  OF  VARIANCE  FOR  LASA  MAGNITUDES 


Thus  far,  simple  models  have  been  used  to  make  inferences  from  LASA  ampli- 
tudes or  magnitudes  concerning  the  bias  of  LASA  relative  to  NEIS,  the  nature 
of  the  differences  between  the  array  beam  and  subarray  beams,  and  structural 
effects  causing  scatter.  A relatively  large  standard  deviation,  0.15  m^ 
existed  for  magnitudes  recorded  over  the  13  LASA  subarrays  studied.  In  this 
section,  to  describe  this  scatter  and  to  estimate  expected  amplitude  variation 
among  nearby  receivers,  a more  complex  model  will  now  be  fit  to  this  LASA 
data.  The  results  can  be  applied  to  predictions  of  magnitude  variation  at  a 
single  station  using  closely-spaced  sources. 

Fixed  Effects  Model  for  LASA  Magnitudes 

Assume  the  following  model  to  describe  the  LASA  magnitude  data: 


LASA  NEIS 
m.  - m..  = p 

ijk  jk 
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where 
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= magnitude  measured  at  the  i-th  subarray  for  the  k-th  event 
in  sector  j. 
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= magnitude  reported  by  NEIS  for  the  k-th  event  in  sector  j, 

= mean  LASA  magnitude  bias  relative  to  NEIS  magnitude  for  all 
events, 

= i-th  subarray  effect, 

= j-th  sector  (azimuth)  effect, 

= interactions  between  i-th  subarray  and  j-th  sector, 

= k-th  event  effect  within  sector  j, 

= a random,  normally-distributed  noise  term  with  zero  mean  and 
variance  o^. 


e 

This  is  a fixed-effect  model  with  interactions;  a detailed  presentation 
of  the  model  can  be  found  in  Scheffe  (1959).  Note  that  6^  terms  are  not 
interactions,  but  rather  nested  effects  within  each  j-th  sector. 


Scheffe,  H. , 1959.  The  Analysis  of  Variance.  New  York,  John  Wiley  and  Sons, 
Inc. 


Although  equalizing  the  number  of  events  from  each  sector  was  not  necessary 
to  arrive  at  an  analytical  solution  for  the  various  terms,  we  used  only  events 
in  the  data  base  recorded  properly  at  all  subarrays  for  each  sector.  To 
reduce  the  computations,  the  number  of  events  in  a sector  was  kept  under  fifty. 
This  restriction  resulted  in  no  significant  loss  of  precision  in  the  estimated 
terms,  and  most  sectors  had  initially  fewer  than  fifty  events  (Appendix  A). 

For  each  j-th  sector,  then,  the  k index  runs  from  1 to  a variable  k^.  To 
arrive  at  independent  solutions  for  the  terms,  these  essential  side  conditions 
were  imposed: 

a =e  =6.  m y . = =0 


where  the  dot  notation  means  averaging  over  the  index  it  replaces.  With 
these  conditions,  the  least-square  estimates  are  obtained  as 

A 

u = y 

• • • 

= y±  - y 

sj  ‘ y.i.  ‘ y... 

in  - y«.  - yi..  - y.j. +y... 

5jk  y.jk  y.j. 


The  estimates  for  4 and  the  a^'s,  g^'s,  and  A„'s  are  listed  in  Table  V. 
Based  upon  initial  estimates  of  the  above  parameters  in  the  model,  the  e 


term 
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eijk  ’ yijk  - “i  - yu  - sj  - 5Jk 


ijk 
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was  used  to  identify  large  measurement  errors  in  the  data  base.  Although 

A 

each  observation  with  l^jjJ  > 0.3  was  rechecked,  only  a few  readings  were 
found  in  error,  and  these  were  eliminated  or  changed  before  computing  final 
estimates  of  the  model  parameters. 

According  to  Scheffe's  (1959)  work,  an  analysis-of-variance  framework 
can  be  constructed  for  the  particular  model  shown  in  Table  VI.  Tests  for 
significance  of  the  various  terms  in  the  model  are  in  the  form 


SS  /DF 
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model  for  LASA  magnitudes. 

SOURCE  OF  VARIANCE  SUM  OF  SQUARES  DEGREES  OF  FREEDOM 


of  sectors  (10) 


where  x stands  for  a,  6,  6,  or  y and  DF^  for  the  degrees  of  freedom  associated 
with  parameter  x.  and  Nq  are  the  number  of  earthquakes  and  observations, 

respectively,  thus,  Nq  * x 13  subarrays.  The  F-statistic  test  was  exceeded 
for  each  of  the  terms  at  the  a = .99  level  of  confidence,  making  that  model 
appropriate  for  LASA  magnitudes.  Contrasts  among  main  effects  for  the  sub- 
arrays can  be  tested,  an  exercise  showing  that  the  values  for  proximate 
subarrays  may  be  closer  in  value,  but  not  significantly  different,  than  for 
those  of  more  distant  subarrays.  Using  Scheffe’s  T-method  of  comparison,  con- 
trasts among  the  subarray  effects  in  the  D,  C,  and  even  B rings  were  estab- 
lished to  be,  at  a high  confidence  level,  non-zero  in  most  cases,  even  though 
the  range  in  the  ct^'s  is  small.  Even  greater  confidence  exists  for  contrasts 
among  sector  effects  Bj  because  of  their  larger  range. 

Event  effects  5 ^ were  not  listed  in  Table  V,  but  they  were  generally 
large,  ranging  from  nearly  -1  to  +1  magnitude  unit.  Assuming  that  NEIS  mag- 
nitudes are  unbiased  estimates,  interpretations  of  event  effects  6 ^ must,  in 
terms  of  source  radiation  pattern  and  geophysical  anomalies,  be  peculiar  to  the 
exact  source-LASA  path.  Although  subdividing  the  10  sectors  into  more  numer- 
ous event  regions  would  reduce  the  scatter  of  these  terms  within  a sector, 
this  refinement  was  unwarranted  because  of  an  increase  in  degrees  of  freedom. 


The  p term  at  LASA  has  an  overall  -0.09  magnitude  bias  for  seismic  events 
within  100°  around  the  globe,  a figure  in  good  agreement  with  this  study's 
earlier  analysis.  The  term  represents  the  mean  bias  of  subarray  magnitudes 
relative  to  p.  Their  small  range  implies  that  a gross  station  effect  (or  cor- 
rection term)  will  vary  only  slightly  over  the  dimension  of  LASA,  roughly  50 
km,  and  that  station  effects  are  not  the  cause  of  the  large  magnitude  scatter 
at  LASA.  The  larger  3.  effects  are  azimuthally-dependent  gross  correction 
terms  due  to  the  structure  in  the  crust  and  upper  mantle  under  LASA.  However, 
a portion  of  these  effects  can  be  attributed  to  consistent  source  or  path 
bias  relative  to  LASA  for  events  within  the  delineated  sectors.  For  example. 
Sector  7 is  dominated  by  earthquakes  in  the  Fij i-Tonga-Kermadec  region,  where 
the  fault-plane  strikes  have  fairly  consistent  azimuths  relative  to  LASA. 

The  interactions  in  Table  V are  random  and  have  a range  of  nearly 

0.5  m^  unit.  These  terms  are  closely  related  to  the  effects  on  signal  wave 
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forms  from  the  laterally  inhomogenous  structure  beneath  LASA  and  they  absorb 
much  of  the  scatter  represented  by  the  high  values  resulting  when  equation 
(8)  was  applied. 


The  noise  term  in  our  model,  eijk»  *s  statistically  represented  by  its 

estimated  variance,  computed  as  .011  using  the  form  in  Table  V;  thus  a = 

0.104.  Assuming  a normal  distribution  for  the  noise  term  implies  that  for 

5%  of  the  observations  the  error  exceeds  an  absolute  value  of  0.21  m,  unit. 

b 

Thus,  even  a detailed  model  does  not  adequately  represent  observed  magnitude 
at  a LASA  subarray,  if  a confidence  factor  of  two  in  the  predicted  value  is 
required.  Further,  this  result  ignores  the  large  and  capricious  values  of 
the  6jk's,  which  are  mostly  unpredictable  and  add  as  much  as  ±1.0  to  the 
uncertainty  in  absolute  source  magnitude. 

Nature  of  Station  Magnitude  Corrections 

Several  investigators  have  already  obtained  seismic  station  magnitude 
corrections  independent  of  azimuth  (Cleary,  1967;  Carpenter  et  al.,  1967; 
Evernden  and  Clark,  1970;  Booth  et  al.,  1975;  North,  1977).  The  results  pre- 
sented here  for  LASA  provide  some  gauge  of  the  utility  of  such  corrections 
when  applied  to  any  given  event.  Consider  each  LASA  subarray  as  a distinct 
station  with  a correction  given  by 

Cij  = U + Oi  + Y±j  + Bj#  (12) 

This  definition  includes  the  sector  effect  and  its  interaction  with  each  sub- 
array and  it  is  a more  precise  station  correction  than  a gross  single-dimen- 
sion station  term.  Although  subdividing  the  sectors  into  small  source  regions. 


Cleary,  J.,  1976.  Analysis  of  the  amplitude  of  short-period  P-waves  recorded 
by  Long  Range  Seismic  Measurements  stations  in  the  distance  range  30  to 
102  degrees,  J.  Geophys.  Res.,  72,  4705-4712. 

Carpenter,  E.  W.,  P.  D.  Marshall,  and  A.  Douglas,  1967.  The  amplitude-dis- 
tance curve  for  short-period  teleseismic  P-waves,  Geophys.  J.  R.  Astr. 
Soc.,  13,  61-70. 


Evernden,  J.  F.,  and  D.  M.  Clark,  1970.  Study  of  teleseismic  P.  II-Amplitude 
data,  Phys.  E.  Planet.  Int.,  4^,  24-34. 

-35- 


F 


as  Berteussen  and  Husebye  (1974)  did  with  NORSAR  amplitude,  can  achieve  fur- 
ther precision,  the  present  limited  data  base  will  not  allow  such  refinement 
for  most  sectors.  Recall  that  the  main  effects  and  B.»  and  their  inter- 
actions, were  found,  at  very  high  confidence  levels,  to  be  significantly  dif- 
ferent from  zero.  Thus,  plots  of  C„  in  Figure  6 for  all  the  subarrays  can 
be  presented  as  meaningful  data.  In  each  case,  connected  sector  points  have 
a positive  or  negative  residual  relative  to  the  NEIS  reference  set  represented 
by  the  dashed  circle,  with  radius  equal  to  one  magnitude  unit,  and  an  angle 
equivalent  to  the  center  azimuth  of  each  sector.  In  this  figure  important 
phenomena  represented  include  the  large  azimuthal  variations  at  any  given  sub- 
array and  the  imperfect  correlation  of  patterns  among  the  subarrays  that  span 
an  area  roughly  50  km  in  diameter.  Evidently,  "station  correction"  is  a 
highly  variable  quantity. 

If  from  equation  (9)  a prediction  of  event  magnitude  from  a single  obser- 


vation is  given  by 


NEIS . , LASA 

m.  . (l)  = m. 
ij  ijk 


- p -a.  - y . . - 6 

1 iJ  j. 


then  the  difference  between  prediction  and  actual  observation  is 

NEIS  NEIS..  J 

mjk  " mjk  = 6jk  + eijk  M 

Note  that  6 ^ cannot  be  known  a_  priori  and,  therefore,  cannot  be  used  to  cor- 
rect station  magnitudes  before  a mean  event  magnitude  is  estimated  from  many 
stations.  From  this  analysis  of  LASA  magnitudes,  6 has  a range  of  roughly 

3 k. 

two  magnitude  units  and,  e, . , with  a standard  deviation  of  approximately  0.1, 

K * 

is  negligible  in  comparison.  Thus,  , which  Figure  6 shows  has  a spread  of 
less  than  one  magnitude  unit,  accounts  for  no  more  than  one-half  the  observed 
magnitude  difference  between  a LASA  subarray  prediction  and  the  NEIS  mag- 
nitude. This  result  implies  that,  although  using  azimuthally-varying  correct- 
ions can  produce  significant  gains  in  precision  of  a stations’s  magnitude 
estimation,  a significant  portion  of  unpredictable  error  will  still  exist 
unless  more  detailed  or  regionalized  corrections  are  developed. 


Berteussen,  K.  A.,  and  E.  S.  Husebye,  1974.  Amplitude  pattern  effects  on 
NORSAR  P-wave  detectability,  NORSAR  Scientific  Report  No.  1-74/75, 
Kjeller,  NORWAY. 
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Figure  6.  Azimuthal  variations  of  LASA  subarray  magnitude  corrections.  Dash- 
ed circle  in  each  subarray  represents  zero  mb  correction.  Inside 
the  circle  smaller  subarray  m^  is  represented;  thus,  positive  correc 
tions  should  be  applied,  and  vice-versa  for  points  outside. 


Precision  of  a Single  Site  in  Estimating  Source  Magnitude 

Although  predicting  seismic  magnitude  at  one  site  from  a nearby  site  is 
of  some  interest,  the  pattern  of  LASA  magnitudes  can  also  be  used  to  investi- 
gate the  more  practical  need  to  predict  the  magnitude  (or  its  uncertainty)  at 
a single  site  when  the  source  position  is  varied.  This  situation  arises  in 
any  "calibration"  context  where  information  gained  from  one  seismic  event  is 
used  to  correct  the  observations  of  one  nearby.  In  such  cases,  both  accuracy 
of  the  calibration  and  the  level  of  confidence  of  the  predictions  are  impor- 
tant. 

If  seismic  reciprocity  is  invoked  (Knopoff  and  Gangi,  1959),  analyzed 
LASA  magnitude  data  is  useful  in  approximating  anticipated  amplitude  variations 
either  at  a single  receiver  site,  when  the  source  is  moved  over  an  area  roughly 
the  extent  of  the  LASA  array,  or  on  the  order  of  five  to  ten  wave-lengths  for 
signals  in  the  routine  LASA  detection  passband  of  0. 9-1.4  Hz.  Let  m^^  and 
be  two  subarray  magnitudes,  consistent  with  the  definition  of  equation 
(9).  At  any  subarray,  these  magnitudes  over  the  event  set  have  no  easily 
formulated  statistical  distribution.  However,  if  we  consider  the  normalized 
form 

-r  A /-l  A VTT?  TP  A A A A * 

(15) 
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J k j l ij  k 


by  this  equation,  the  left-hand  side  can  be  approximated  by  a variable  with 
a normal  distribution  having  mean  ou  and  a variance  given  by  the  randomness 
Y^j  + e^jjc  over  all  events.  We  now  adopt  the  simplified  notation  of 


LASA  NEIS  * * 

m.  - m.  — m..  — o,t  — p.  — u 

l ijk  jk  jk  pj  M 


(16) 


for  a subarray's  normalized  magnitude  and  then  express  the  conditional  distri- 
bution of  m",  for  a suite  of  events,  as 


f(m  , m ) = ■ ■ o 

m Tar  n ~2 


exp 


IV 


2(l-p.^)  o^ 


(17) 


Knopoff,  L. , and  A.  F.  Gangi,  1959.  Seismic  Reciprocity  Geophysics,  24  (4) 
681-691. 
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where  the  mean  $ is  given  by 

^ii'  = pii"  (mi  “ “±)  + a±'  (18) 

Equation  (17)  is  the  relation  appropriate  to  jointly-distributed,  normal,  cor- 
related variables.  To  determine  we  compute  the  covariance  of  nu  and 


(m 


LASA 

ijk 


NEIS 

m., 

3k 


LASA  NEIS 

i. - m.. 
i jk  Jk 


(19) 


and  the  variance  of  m 


i 


as 


and  similarly  for  0^2 


_1 

K. 

3 


LASA  NEIS  . * * 

’I  (mijk-mjk-  6jk'  Bj-y"mi..)2  <20> 

The  squared  correlation  coefficient  is  defined  as 


°i  V (21) 

Aki  (1973)  has  already  shown  that  p^,  is  dependent  upon  the  subarray  spacing 
at  LASA  for  spectral  amplitudes  near  1 Hz,  a finding  consistent  with  Chernov's 
(1960)  theory.  Figure  7 shows  the  p ^ computed  for  the  13  subarrays  of  this 
study.  The  trend  is  toward  lesser  correlation  with  increasing  separation, 
which  is  similar  to  Aki's  results  because  magnitudes  measured  in  this  study 
were  associated  with  periods  predominantly  near  one  second.  The  o^/l-p.^2, 
values  are  also  plotted  which  are  the  standard  deviation  of  the  conditional 
distribution  of  magnitude  at  one  subarray,  given  that  at  another.  These 
values  are  in  part  reflections  of  results  displayed  in  Figure  6 on  the  ampli- 
tude patterns  of  the  LASA  subarrays. 


How  the  principle  of  seismic  reciprocity  between  sources  and  receivers 
(Knopoff  and  Gangi,  1959)  is  related  to  equality  of  amplitude  variances  is 
discussed  in  Appendix  B.  If  reciprocity  holds,  then  the  correlations  and 
uncertainties  expressed  in  Figure  7 can  be  converted  into  appropriate  uncer- 
tainties for  an  array  of  sources  rather  than  receivers.  For  a pair  of  iden- 
tical sources  within  ten  km  of  each  other,  roughly  the  smallest  inter-subarray 
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Figure  7.  Estimates  of  signal  amplitude  correlations,  variances,  and  differ- 
ences versus  receiver  separation  using  LASA  subarrays.  By  seismic 
reciprocity,  this  is  an  estimate  of  variations  to  be  expected  at  a 
single  site  due  to  similar  changes  in  source  position. 


distance  at  LASA  (between  AO  and  B ring),  correlation  of  the  amplitudes 
received  at  a single  fixed  site  should  be  near  0.5.  Corresponding  to  this 
correlation  the  standard  deviation  of  the  measured  magnitude  for  the  second 
source  minus  that  of  the  first  source  should  be  near  0.10.  For  maximum 
separation  of  sources  corresponding  to  the  LASA  D-ring  spacing,  the  correla- 
tion falls  to  around  zero,  and  the  standard  deviation  of  the  measured  magni- 
tude of  the  second  source,  relative  to  that  of  the  first  one,  would  be  near 
0.15.  Note  that  the  variation  in  the  a.  is  not  considered  here  because  it 

l 

incorporates  the  specific  peculiarity  of  the  LASA  structure  and  it  cannot  be 
projected  into  a statistical  consideration  of  source  regions.  Its  values 
(see  Table  V)  were  small  at  LASA  and  including  it  would  add  no  more  than  0.1 
to  the  uncertainty  for  the  largest  separations  and  much  less  for  the  smallest 
separations  considered. 

These  results  imply  that  magnitude  calibration  for  sources  has  a very 
limited  range,  probably  less  than  is  commonly  thought  satisfactory  for  a cal- 
ibration of  travel-time  residuals,  except  in  dipping-plate  source  regions. 

If  LASA  structure  is  similar  to  the  crust  and  upper  mantle  in  earthquake 
source  regions,  in  terms  of  expected  variations  in  magnitude,  then  a magnitude 
calibration  for  a given  source-receiver  pair  cannot  be  expected  to  remain 
within  0.3  m^  unit  with  better  than  95%  confidence  when  a new  source  occurs 
just  50  km  from  the  calibration  source.  Clearly,  an  advantage  exists  in 
multiple  receivers  because  the  improvement  through  network  calibration  should 
behave  as  the  square  root  of  the  number  of  receivers. 

The  question  of  whether  the  conditions  for  seismic  reciprocity  are  valid 
for  the  LASA  array  of  receivers  is  not  studied  in  this  report.  If,  as  Bala- 
chandran  (1974)  suggests  for  a smaller-scale  experiment,  the  LASA  seismometers 
are  within  a few  wavelengths  of  the  inhomogeneities  that  give  rise  to  the  am- 
plitude and  phase  fluctuations  found  in  this  and  numerous  other  studies,  then 
reciprocity  may  not  hold.  The  scatter  of  amplitudes,  enhanced  by  secondary 
wave  fields,  over  the  LASA  array  may  ;»ive  too  large  an  estimate  of  the  scatter 
to  be  expected  in  the  reciprocal  case  of  an  array  of  sources. 


Balachandran,  K. , 1974.  Noninterchangeability  of  sources  and  receivers  in  a 


BEHAVIOR  OF  LASA  AMPLITUDES  AS  A FUNCTION 
OF  SENSOR  AND  SOURCE  SEPARATION 

This  section  focusses  on  variations  in  receiver  amplitudes  as  a function 
of  increasing  source  separation.  The  discussion  clarifies  nearly  all  the 
important  results  of  this  report  that  pertain  to  statistical  fluctuations  of 
LASA  amplitudes. 

Consider  the  log  amplitudes  a^k  at  two  subarrays  i and  1'  for  events  k 
and  k , with  k indexing  over  all  events  and  not  just  those  in  a particular 
sector.  Then  the  quantity 

Aa  = (aik  - a^)  - (a^  - Vk*> 


measures  the  imprecision  of  the  calibration  of  the  subarray  differences  from 
one  event  to  another.  As  already  discussed,  this  quantity  should  be  an  increas- 
ing function  of  subarray  separation  and  of  event  separation.  Here  we  compute 
the  standard  deviation  of  Aa  versus  increases  in  both  these  separations  thus. 


o[Aa( jr ) ,dy)]= 


N faik  aj"k)  ~ (a 
N-l 


ik' 


ai'k'^ 


(22) 


where  the  summations  is  made  over  all  subarray-event  combinations  falling 
within  given  incremental  values  of  subarray  separation  |r|  and  central  angle 
separation  dy.  The  dy  is  calculated  as 

dy  = arccos  (X^X^  + + W') 

where 


*k  = Sin  ik  * 
Yk  = Sln  *k  * 
Zk  = COS  ik 


cos  0, 
k 

sin  0, 
k 


and  similarly  for  event  k'.  The  angles  ik  and  0k  are  the  incidence  angle  and 
back  azimuth  of  the  signal  arriving  at  LASA  from  event  k.  So  dy  is  a measure 
of  the  central  angle  between  the  two  rays  that  impinge  on  an  imaginary  hemi- 
sphere surrounding  LASA.  The  epicentral  distance  and  depth  of  the  events, 
listed  in  Appendix  I,  are  converted  to  an  incidence  angle  by  linear  inter- 
polation of  the  tables  in  Pho  and  Behe  (1972). 
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When  computing  standard  deviations  according  to  equation  ( 22 ),  subarray 
separations  were  pooled  in  10  km  steps  and  dy  pooled  into  bins  of  < 1°,  < 2°, 

< 4°,  etc.  The  truncated  data  base  defined  in  the  previous  section  was  used 
as  a means  to  avoid  biasing  the  results  with  large  concentrations  of  data 
from  a few  highly  seismic  areas.  The  epicenters  were  further  limited  to  A > 
30°  from  LASA  because  the  i^  vs.  A relation  is  highly  dependent  upon  crustal- 
upper  mantle  structure  at  short  distances.  Also,  epicenters  with  A < ^0° 
were  discarded  because  highly  variable  fluctuations  were  associated  with  even 
close  epicenters  because  of  core  diffraction  and  plate  effects  on  many  of  the 
events  in  this  range  from  LASA  (see  Figure  1).  With  the  167  remaining  events, 
a(Aa)'s  shown  in  Figure  8 were  computed.  Each  point  represen  s no  fewer  than 
N=500  combinations  summed  in  equation  (22).  The  pattern  here  is  relatively 
smooth,  increasing  from  roughly  a = .10  for  small  |r|  and  small  dy  (in  gen- 
eral agreement  with  results  in  previous  sections)  to  nearly  o = .40  for  the 
largest  values  of  both.  Therefore,  results,  shown  in  Figure  8 suggest  that 
magnitude  calibration  is  not  precise  for  neighboring  events  and  that  it 
becomes  almost  meaningless  for  events  distant  from  one  another.  Note  that  a 
given  dy  maps  into  values  of  inter-event  distance  which  are  dependent  upon  the 
relation  of  the  event  locations  to  the  great-circle  path  back  to  the  receiver. 
For  example  at  A = 90°,  a 1°  difference  in  dy  translates  to  a 1°  (*•  110  km) 
difference  in  epicenters  perpendicular  to  the  path  but  to  roughly  a 10° 

(~  1100  km)  difference  in  epicenters  parallel  to  the  path.  Thus,  magnitude 
calibration  could  be  applicable  to  much  larger  source  location  changes  if 
they  are  along  the  raypath  direction  rather  than  normal  to  it.  However,  cal- 
culations needed  to  precisely  verify  this  principle  were  not  performed. 

Results  similar  to  those  for  LASA  in  Figure  8 have  been  obtained  for 
NORSAR  using  the  data  in  Berteussen  and  Husebye  (1974),  but  they  are  not 
presented  here. 


Pho,  T.-T. , and  L.  Behe,  1972.  Extended  distances  and  angles  of  incidence 
of  P-waves,  Bull.  Seism.  Soc.  Am.,  62 , 885-902. 
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CONCLUSION 

LASA  magnitude  bias  relative  to  NEIS,  when  averaged  in  all  directions 
for  the  whole  array,  as  well  as  individual  subarrays,  is  near  -0.1  m^.  How- 
ever, when  viewed  in  azimuthal  sectors,  LASA  m^  biases  reveal  large  coherent 
variations,  indicating  that  generally  accepted,  single-value,  station  bias 
terms  cannot  reduce  bias  to  negligible  proportions  for  single  stations  mag- 
nitude observations,  even  if  they  are  determined  to  be  statistically  signi- 
ficant with  a large  data  base.  Subarray  azimuthal  biases  show  uniform  vari- 
ation among  subarrays,  suggesting  that  the  effect  of  regional  bias,  possibly 
related  to  upper-mantle  attenuation,  indeed  exists  at  LASA.  Interpretations 
of  regional  bias  suggests  that  LASA  is  located  at  the  edge  of  the  anomalous 
Western  United  States  structure.  However,  local  crustal  and  upper-mantle 
heterogeneity  apparently  has  an  even  greater  effect  than  regional  structure 
on  the  m^  bias  of  each  subarray.  The  result  is  that  local  m^  biases  also 
reveal  characteristic  features  of  crustal  and  upper-mantle  structure  and 
they  show  some  linear  correlation  with  values  of  travel-time  residuals. 

From  the  linear  regression  analysis  standard  deviations  of  subanay  mag- 
nitudes were  found  to  be  at  the  0.15  m^  level  within  each  sector.  Using  the 
fixed  effects  model,  these  equivalent  values  are  represented  in  the  sector- 
subarray  interaction  term  > which  also  showed  a similar  range  of  fluctua- 
tions. Further,  the  event  effect  term  6 ^ gives  as  much  as  much  as  1.0  m^ 
uncertainty  in  source  magnitude.  Subarray  effects  were  small  when  averaged 
over  all  events,  and  so  they  are  not  the  cause  of  the  large  LASA  magnitude 
scatter. 

Within  the  context  of  station  calibration,  the  confidence  level  of  pre- 
dicting the  m^  at  a station  for  an  event  from  calibrations  based  on  another 
event  occurring  in  a nearby  source  area  was  examined.  Assuming  seismic  reci- 
procity results  showed  that  the  correlation  of  a single  receiver's  magnitudes 
for  repeated  observations  of  two  events  occurring  only  50  km  apart  from  each 
other  would  be  near  zero.  This  result  implies  that  magnitude  calibration  for 
sources  is  very  limited  in  space,  probably  to  an  area  much  less  than  that 
believed  to  hold  for  the  calibration  of  travel-time  anomalies.  Uncertainty 
is  a factor  of  two  (±0.3  m^) , at  the  95%  confidence  level,  in  predicted 
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amplitude  at  a station  for  a source  moved  only  50  km  from  the  calibration 
source  in  a geologic  area  of  complexity  similar  to  LASA's.  For  sources 
several  tens  of  degrees  apart  the  magnitude  uncertainty  at  the  95%  confidence 
level  is  +0.8  m^  for  stations  50  km  apart  in  regions  similar  to  LASA.  To  the 
extent  that  the  crustal  effects  discussed  in  this  report  may  be  regarded  as 
"random"  in  a network  of  stations,  they  may  be  miminized  by  network  averaging. 
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APPENDIX  A 

List  of  events  and  event  parameters  used 
in  this  study,  grouped  in  sectors  L01  to  L10. 


LA5A  MAGNITUDE,  DMB=MLASA-HNEIS, SECTOR  L01 


3? 

EVFNTS 

ARE  ANALYZED  FROM 

350.1 

TO  25.0 

AZI 

EPX  NO. 

DATE 

MNEIS 

MLASA 

HL-HN 

DIST 

AZI 

DEPTH 

BFACTOR 

11680 

740404 

5.  10 

4.98 

-0.  12 

94.  1 

1.0 

21.0 

3.80 

17200 

740406 

5.20 

5.44 

0.24 

96.2 

1.0 

84.0 

4.00 

44000 

740122 

5.00 

5.20 

0.20 

93.2 

1.0 

50.0 

3.80 

71405 

740219 

5.00 

4.70 

-0.  30 

97.3 

2.0 

05.0 

4.00 

7252^ 

740220 

5.  10 

5.44 

0.34 

92.6 

0.1 

25.0 

3.90 

46930 

740426 

5.00 

4.58 

-0.42 

97.3 

2.0 

80.0 

4.00 

73035 

740813 

5.50 

5.65 

0.15 

96.8 

2.0 

208.0 

3.80 

80906 

740517 

5.30 

5.64 

0.34 

96.8 

2.0 

208.0 

3.80 

4180 

740603 

5.  30 

5.76 

0.46 

96.4 

1.0 

100.0 

4.00 

4590 

740603 

5.00 

5.10 

0.10 

94.  1 

1.0 

60.0 

3.00 

8785 

740606 

5.30 

5.43 

0.13 

96.9 

2.0 

214.0 

3.80 
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740705 

5.10 

5.09 

-0.01 

95.9 

1.0 

33.0 

4.00 
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740723 

5.00 
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33.0 

3.80 
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5.30 

4.51 
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97.7 

1.0 
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33.0 
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0.17 
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0. 1 

33.0 

3.80 
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96.7 

1.0 
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3.00 

s 3905 

740929 

5.  20 

5.15 

-0.05 

65.7 

4.0 

0.0 

3.70 

13815 

741003 

5.00 

4.93 

-0.07 

06.8 

2.0 

198.0 

3.80 

3985 

741113 

5.10 

5.02 

-0.00 

87.4 

19.0 

42.0 

3.70 

1 2 305 

741117 

5.20 

5.48 

0.28 

98.7 

15.0 

43.0 

4.20 

28345 

741127 

5.00 

5.47 

0.47 

94.2 

22.0 

50.0 

3.85 

87460 

741230 

5.30 

5.12 

-0.18 

97.2 

3.0 

116.0 

4.00 

51926 

740130 

5.40 

5.  78 

n.  38 

03.4 

357.0 

0.0 

3.70 

61440 

740710 

5.30 

5.  39 

0.09 

83.4 

357.0 

0.0 

3.70 

17285 

7408 1 1 

5.10 

5.14 

0.04 

94.0 

359.0 

12.0 

3.80 

17475 

740811 

5.10 

4.91 

-0.19 

94.  1 

359.0 

29.0 

3.80 

20170 

740812 

5.20 

4.75 

-0.41 

94.1 

359.0 

27.0 

3.80 

50300 

740827 

5.30 

5.40 

0.10 

93.9 

359.0 

, 33.0 

3.80 

77165 

740912 

5.20 

5.13 

-0.07 

94.  1 

359.0 

33.0 

3.80 

7?285 

740913 

5.20 

5.40 

0.20 

83.4 

357.0 

0.0 

3.70 

7768 

740929 

5.40 

5.97 

0.57 

92.8 

356.0 

33.0 

3.90 

42750 

741016 

5.50 

5.76 

0.26 

83.2 

356.0 

0.0 

3.70 

8650C 

741103 

5.20 

5.  39 

0.19 

89.4 

354.0 

33.0 

3.70 

A-l 


jBMKMttWWttS8® 


15  BIST  QUALITY  JRACTlC^a 


xo  ddC- 


LASA  MAGNITUDE,  DM P=MLA SA-HNEIS , S ECTOR  L02 


22 

EVENTS 

ARE  ANALYZED  FROM 

25.1 

TO  65.0 

AZI 

EPX  NO. 

DATE 

HNEIS 

MI.ASA 

ML-NN 

DIST 

AZI 

DEPTH 

BFACTOR 

85505 

741229 

5.20 

5.06 

-0.14 

48.  4 

35.0 

33.0 

3.60 

6«9R0 

741214 

5.00 

4.48 

-0.52 

82.3 

39.0 

11.0 

3.60 

63290 

741214 

5.  30 

5.  96 

0.66 

82.7 

39.0 

33.0 

3.70 

10705 

741116 

5.00 

5.32 

0.32 

46.2 

53.0 

33.0 

3.50 

6380 

741114 

5.00 

*.11 

0.11 

83.5 

37.0 

24.0 

3.70 

6350 

741114 

5.10 

5.22 

0.12 

83.5 

37.0 

3.0 

3.70 

6305 

741114 

5.00 

5.05 

0.05 

83.5 

37.0 

19.0 

3.70 

75625 

741029 

5.10 

5.70 

0.60 

76.5 

37.0 

33.0 

3.60 

42541 

741016 

5.00 

5.27 

0.27 

46.  1 

53.0 

33.0 

3.50 

86205 

740917 

5.20 

4.36 

-0.82 

81.0 

38.0 

17.0 

3.50 

64390 

740904 

5.  10 

5.34 

0.24 

83.6 

46.0 

17.0 

3.70 

73180 

740717 

5.10 

5.78 

0.68 

78.7 

31.0 

145.0 

3.60 

60195 

740709 

5.00 

4.79 

-0.21 

87.3 

34.0 

69.0 

3.70 

47835 

740628 

5.00 

5.41 

0.41 

76.6 

50.0 

33.0 

3.60 

44760 

740625 

5.10 

5.09 

-0.01 

48.  3 

35.0 

33.0 

3.60 

39165 

740622 

5.10 

5.41 

0.31 

81.2 

36.0 

36.0 

3.50 

33550 

740620 

5.10 

4.76 

-0.34 

76.4 

37.0 

33.0 

3.60 

17770 

740612 

5.50 

6.13 

0.63 

46.9 

35.0 

1 3.0 

3.60 

81016 

740517 

5.00 

5.60 

0.60 

46.9 

35.0 

33.0 

3.60 

98820 

740323 

5.00 

4.03 

-0.07 

43.9 

53.0 

33.0 

3.20 

98340 

740322 

5.00 

4.80 

-0.20 

47.1 

26. C 

22.0 

3.60 

53670 

740201 

5.20 

5.66 

0.46 

85.1 

34.0 

29.0 

3.70 

_ — ~ '4UAXTX1Z  MuiWllUUUJ 

oopx  ntmisHi)  i©  doc- 


las  a HAONITODE,  D"B  = HLA5A-f!NEIS,  SECTOR  L03 


16 

EVENTS 

ARE  ANALYZED  FRON 

65.1 

TO  110.0 

AZI 

FPX  NO. 

DATE 

"N  FIS 

N LAS  A 

NL-HN 

DIST 

AZI 

DEPTH 

BFACTOR 

64460 

741210 

5.00 

5.11 

0.11 

51.3 

94.0 

33.0 

3.40 

18425 

741120 

5.00 

5.26 

0.26 

84.  1 

100.0 

33.0 

3.70 

98070 

740924 

5.10 

5.12 

0.02 

60.6 

101.0 

33.0 

3.60 

72780 

740909 

5.00 

5. 10 

0.  10 

52.4 

76.0 

33.0 

3.40 

19575 

740812 

5.00 

5.04 

0.04 

89.8 

90.0 

33.0 

3.70 

50530 

740701 

5.00 

4. Op 

-0.02 

78.9 

99.0 

33.0 

3.50 

47391 

740627 

5.30 

5.39 

0.09 

78.8 

99.0 

33.0 

3.50 

46991 

740627 

5.40 

5.48 

0.08 

78.9 

99.0 

33.0 

3.50 

89755 

740523 

5.10 

4.73 

-0.37 

51.6 

89.0 

33.0 

3.40 

77243 

7405 1 c 

5.00 

4.47 

-0.  53 

51.5 

89.0 

33.0 

3.40 

62975 

740507 

5.20 

4.71 

-0.49 

61.5 

102.0 

33.0 

3.70 

31475 

740417 

5.10 

5.00 

-0.  10 

52.9 

75.0 

33.0 

..  3-“9 

6275 

740330 

5.00 

5.28 

0.28 

76.2 

99.0 

33.0 

3.60 

66898 

740214 

5.40 

5.00 

-0.  40 

55.2 

94.0 

33.0 

3.50 

6165 

740604 

5.00 

5.37 

0.37 

64.2 

102.0 

33.0 

3.70 

4Q0  30 

740127 

5.  10 

4.94 

-0.  16 

51.6 

78.0 

33.0 

3.40 

A- 3 


JBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

ZBQM  Ouri  iutaasHED  to  dlc 


L ASA  MAGNITUDE,  PMB=MLA SA-NNEIS, SECTOR  L04 


6 

EVENTS 

ARE  ANALYZED  FROM 

110.1 

TO  130.0 

AZI 

EPX  MO. 

DATE 

HMEIS 

MLASA 

Hl-HN 

DIST 

AZI 

DEPTH 

BFACTOR 

17740 

741120 

5.20 

5.07 

-0.13 

50.6 

113.0 

33.0 

3.40 

47385 

741018 

5.20 

5.48 

0.  28 

46.3 

113.0 

45.0 

3.50 

24145 

741006 

5.10 

5.07 

-0.03 

46.6 

113.0 

53.0 

3.60 

66640 

740214 

5. 10 

4.89 

-0.21 

40.0 

119.0 

7.0 

3.10 

75895 

741029 

5.10 

5.16 

0.06 

51.1 

120.0 

33.0 

3.40 

40835 

740118 

5.30 

5.02 

-0.28 

41.0 

120.0 

82.0 

3.20 

IBIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
7R0H  OUEY  fUMUSHED  TO  DCC  — 


LASA  MAGNITUDE,  DM B= ML A SA -MNEIS , S ECTOP  L05 


97 

EVENTS  . 

ARE  ANALYZED  FROM 

140.  1 

TO  160.0 

A7.I 

EPX  NO. 

DATE 

MNEIS 

MLASA 

ML-MN 

DTST 

AZI 

DEPTH 

BFACTOR 

66  060 

741025 

5.50 

6.38 

0.88 

32.7 

156.0 

120.0 

3.40 

95060 

740921 

5.30 

6.07 

0.77 

77.8 

144.0 

106.0 

3.60 

20746 

741007 

5.00 

4.87 

-0.  13 

64.  1 

149.0 

33.0 

3.70 

92655 

740920 

5.00 

4.70 

-0.  30 

35.3 

153.0 

33.0 

3.40 

15960 

74  061  1 

5.10 

5.43 

0.33 

82.7 

146.0 

132.0 

3.70 

86021 

740726 

5.00 

5.  43 

0.43 

70.5 

143.0 

157.0 

3.60 

20600 

740409 

5.10 

5.36 

0.26 

84.4 

150.0 

94.0 

3.70 

47315 

741206 

5.40 

5.61 

0.21 

43.3 

145.0 

46.0 

3.20 

37485 

740920 

6.40 

5.91 

0.51 

82.3 

147.0 

93.0 

3.60 

28785 

740816 

5.00 

5.3  6. 

0.36 

43.2 

144.0 

36.0 

3.20 

90705 

741231 

5.40 

5.65 

0.25 

34.7 

155.0 

75.0 

3.40 

86605 

741229 

5.50 

5.82 

C.  32 

86. 1 

150.0 

99.0 

3.60 

70720 

741217 

5.00 

5.29 

0.29 

36.7 

150.0 

50.0 

3.20 

62540 

741213 

5.30 

5.66 

0.36 

73.5 

143.0 

118.0 

3.50 

35405 

740819 

5.20 

6.07 

0.87 

37.3 

151.0 

67.0 

3.20 

42120 

741204 

5.10 

5.64 

0.54 

76.5 

145.0 

53.0 

3.60 

3831  5 

741202 

5.00 

5.64 

0.64 

82.4 

147.0 

92.0 

3.60 

25300 

741124 

5.00 

5.25 

0.  25 

45.  9 

146.0 

36.0 

3.50 

6560 

741114 

c .40 

6.49 

0.0° 

65.0 

148.0 

33.0 

3.70 

2085 

741112 

5.50 

6.08 

0.58 

86.0 

150.0 

90.0 

3.60 

1615 

741112 

5.00 

u.55 

-0.45 

62.3 

149.0 

33.0 

3.70 

93520 

741107 

5.40 

6.25 

0.85 

69.9 

143.0 

153.0 

3.60 

81220 

741031 

5.20 

5.51 

0.31 

69.5 

143.0 

83.0 

3.60 

70  50  5 

741026 

5.30 

5.60 

0.  30 

68.3 

145.0 

110.0 

3.70 

38  1 c 5 

741014 

5.10 

5.14 

0.04 

c5.6 

143.0 

162.0 

3.50 

32810 

741011 

5. 10 

5.  36 

0.26 

83.6 

150.0 

65.0 

3.70 

31180 

741011 

5.30 

5.48 

0. 18 

67.6 

146.0 

99.0 

3.70 

30135 

741010 

5.  30 

5.34 

0.04 

64.4 

148.0 

27.0 

3.70 

14120 

741003 

5.00 

5.04 

0.04 

64.5 

149.0 

33.0 

3.70 

74879 

740911 

5.20 

5.65 

0.  45 

73.5 

143.0 

113.0 

3.50 

60770 

740902 

5.40 

5.87 

0.47 

79.3 

144.0 

163.0 

3.50 

50040 

740827 

5.50 

6.07 

0.57 

52.7 

145.0 

147.0 

3.70 

49055 

740827 

5.40 

4.72 

-0.68 

89.9 

154.0 

23.0 

3.70 

4791  5 

"’40826 

5.10 

5.01 

-0.09 

33.2 

152.0 

3 3.0 

3.40 

45450 

740824 

6.40 

5.24 

-0.16 

85.6 

147.0 

32.0 

3.60 

440«6 

740824 

6.40 

6.81 

0.41 

77.1 

144.0 

109.0 

3.60 

6510 

740^28 

6,00 

6.56 

0.56 

46.2 

146.0 

33.0 

3.50 

63415 

740507 

5.00 

5.  1° 

0.19 

43. 1 

144.0 

3 3.0 

3.20 

80420 

740516 

5.50 

6.58 

0.08 

39.3 

147.0 

36.0 

3.  10 

8105 

74  06  06 

5.10 

5.36 

0.26 

76.  2 

144.0 

116.0 

3.60 

11720 

740600 

5.00 

6.59 

0.59 

77.3 

144.0 

106.0 

3.60 

1 3320 

74  06  09 

5.10 

4.66 

-0.  44 

57.0 

149.0 

52.0 

7.50 

50855 

740701 

5.50 

5.49 

-0.01 

78.3 

141.0 

13.0 

3.60 

62096 

740710 

5.00 

5.24 

0.24 

78.0 

144.0 

121.0 

3.60 

7 3600 

740917 

6.20 

5.50 

0.30 

80.8 

148.0 

3 3.0 

3.50 

76485 

740719 

5.20 

4.96 

-0.24 

85.3 

151.0 

44.0 

3.70 

79920 

740721 

5.00 

4.89 

-0.11 

39.4 

147.0 

78.0 

3.  10 

33335 

74C818 

5.30 

5.44 

0.14 

92.6 

157.0 

21.0 

3.90 

29425 

74  C 1 0 2 

5.30 

6.26 

-0.04 

77.3 

144.0 

99.0 

3.60 

33935 

740107 

5.10 

5.00 

-0.10 

85.6 

151.0 

24.0 

3.60 

39676 

7401 17 

5.30 

5.66 

0.36 

81.4 

148.0 

83.0 

3.50 

43786 

740121 

5.00 

5.26 

0.26 

33.5 

157.0 

33.0 

3.40 

54  080 

740201 

5.20 

5.42 

0.22 

69.4 

147.0 

48.0 

3.70 

55700 

740202 

5.10 

5.40 

0.  30 

34.7 

155.0 

75.0 

3.40 

56515 

740202 

5.20 

5.41 

0.21 

68.8 

147.0 

63.0 

3.70 

66  405 

7U0214 

6.00 

4.98 

-0.02 

91.2 

144.0 

33.0 

3.50 

60811 

740216 

5.30 

5.70 

0.40 

75.5 

145.0 

74.0 

3.60 

72155 

740220 

6.50 

5.84 

0.  34 

93.  9 

147.0 

115.0 

3.70 

73430 

740221 

5.10 

6.81 

0.71 

87.8 

151.0 

87.0 

3.80 

79995 

740228 

6.20 

4.6  3 

-0.57 

41.9 

146.0 

59.0 

3.20 

88720 

740310 

5.10 

4.64 

-0.46 

51.6 

145.0 

43.0 

3.40 

A-5 


IBIS  PAGE  IS  BEST  QUALITY  PBACIICA&LI 

C&r'i  bupaISHED  to  D°Q  ^ — - 


LASA  MAGNITUDE,  DMB= MLA SA -HNEIS , S ECTOR  LOS 

97  EVENTS  ARE  ANALYZED  FROM  140.!  TO  160.0  A7.I 


EPX  NO. 

DATE 

HNEIS 

HLASA 

HL-NN 

DTST 

AZI 

DEPTH 

BFACTOR 

R9482 

740311 

5.  10 

5.26 

0.  16 

61.3 

148.0 

97.0 

3.60 

945 

740324 

5.30 

5.67 

0.37 

86.0 

150.0 

104.0 

3.60 

7605 

740331 

5.00 

6.20 

0.20 

39.  2 

150.0 

49.0 

3.20 

9570 

740402 

5.50 

5.86 

0.36 

85.9 

145.0 

166.0 

3.60 

27466 

7404  14 

5.30 

4.90 

-0.40 

34.  3 

153.0 

138.0 

3.40 

:»it 

740425 

740426 

3:  it 

1:8 

-0.01 

0.13 

71.3 

78.8 

144.0 

«:* 

«:* 

1:18 

7020 

740331 

5.20 

5.39 

0.19 

82.7 

33.0 

3.70 

24355 

92610 

741008 

5.00 

4.72 

-0.28 

55.8 

148.0 

33.0 

3.50 

740823 

5.40 

5.72 

0.32 

79.7 

146.0 

102.0 

3.  40 

4460 

740804 

5.20 

5.47 

0.  27 

73.  4 

146.0 

63.0 

3.60 

45660 

740123 

5.20 

5.73 

0.53 

85.4 

149.0 

115.0 

3.70 

13115 

740107 

5.  30 

5.63 

0.33 

82.7 

149.0 

35.0 

3.70 

80015 

740721 

5.00 

5.20 

0.20 

39.3 

147.0 

97.0 

3.  10 

14880 

740810 

5.20 

5.45 

0.25 

75.5 

143.0 

116.0 

3.60 

79252 

740721 

5.30 

5.34 

0.04 

34.5 

155.0 

70.0 

3.40 

22160 

7404  in 

5.40 

5.31 

-0.09 

34.4 

154.0 

108.0 

3.40 

15  321 

74100a 

5 .00 

5.06 

0.06 

79.3 

140.0 

533.0 

3.50 

70205 

740907 

5.20 

4.67 

-0.  53 

89.8 

155.0 

16.0 

3.70 

21555 

741007 

5.10 

6.07 

-0.03 

39.2 

146.0 

225.0 

3.10 

73820 

741224 

5.40 

6.25 

0.85 

35.  1 

152.0 

155.0 

3.40 

65965 

740905 

5.40 

5.84 

0.44 

49.0 

133.0 

160.0 

3.50 

59895 

740206 

5. 10 

6.44 

0.  34 

48.9 

133.0 

160.0 

3.50 

71435 

7407 1 * 

5.30 

5.28 

-0.02 

46.2 

138.0 

28.0 

3.50 

68250 

740714 

5.20 

5.33 

0.13 

45.9 

138.0 

26.0 

3.50 

67485 

740713 

5.40 

5.38 

-0.02 

46.0 

138.0 

5.0 

3.50 

65206 

740713 

5.50 

5.57 

0.07 

46.3 

138.0 

45.0 

3.50 

65165 

740713 

5.20 

5.25 

0.05 

46.1 

138.0 

23.0 

3.50 

65100 

740713 

c.  30 

5.66 

0.36 

46. 1 

138.0 

24.0 

3.50 

51266 

74  04  2 9 

5.10 

5.41 

0.31 

49.3 

139.0 

87.0 

3. 50 

3776 

740602 

5.2C 

5.53 

0.  33 

48.5 

139.0 

64.0 

3.50 

66100 

740705 

5.20 

4.96 

-0.24 

37.9 

139.0 

23.0 

3.20 

64836 

740713 

5.40 

5.39 

-0.01 

46.  1 

139.0 

33.0 

3.50 

64875 

740713 

5.20 

5.11 

-0.09 

46.6 

139.0 

21.0 

3.60 

65266 

74071  3 

5.10 

4.96 

-0.14 

46.2 

139.0 

33.0 

3.50 

67955 

740713 

5.30 

5.47 

0.17 

46.5 

139.0 

2 3.0 

3.60 

r*  mm p 


IBIS  PJU2S  IS  BBST  QUALITY  PRACtlCA3fi| 
JROM  copy  lyRAiaasL  to  odq 


LAS  A HAGNTT0D5,  n«B=NLA SA-MNFI5, SECTOR  L06 


20 

EVENTS 

APE  ANALYZED  FROM 

160.  1 

TO  210. 

0 AZI 

EPX  NO. 

DATE 

NNFTS 

NLASA 

NL-fIN 

DIST 

AZI 

DEPTH 

BFACTOR 

30830 

740105 

5.00 

4.69 

-0.31 

38.4 

176.0 

33.0 

3.20 

32795 

740107 

6.40 

5.77 

0.37 

49.4 

160.0 

33.0 

3.60 

48275 

740126 

5.10 

4.63 

-0.47 

28.2 

174.0 

33.0 

3.30 

34140 

740418 

5.10 

5.06 

-0.04 

85.7 

170.0 

33.0 

3.60 

34  34  5 

740408 

5.50 

5.58 

0.08 

85.6 

170.0 

33.0 

3.60 

48675 

740428 

5.10 

5.17 

0.07 

50.6 

177.0 

3 3.0 

3.40 

62325 

740*06 

5.00 

4.86 

-0.  14 

77.6 

183.0 

33.0 

3.60 

17500 

740612 

5.00 

5.  16 

0.  16 

30.4 

166.0 

49.0 

3.30 

19640 

74061 3 

5.00 

4.98 

-0.02 

30.7 

164.0 

96.0 

3.20 

43370 

740625 

5.30 

6.09 

-0.  21 

32.5 

160.0 

25.0 

3.40 

55440 

740705 

6.20 

4. 32 

-0.98 

94.3 

16  1.0 

33.0 

3.80 

57010 

740706 

5.10 

5.01 

-0.09 

81.7 

....1.81.0 

13.0 

3.60 

85070 

740725 

5.30 

5.66 

0.  35 

73.7 

187.0 

33.0 

3.50 

20180 

740812 

5. 10 

4.90 

-0.20 

29.5 

168.0 

73.0 

3.30 

59085 

740901 

5.00 

u . 99 

-0.01 

92.  1 

163.0 

33.0 

3.80 

20155 

741006 

5.00 

6.  13 

0.13 

30.  3 

167.0 

51.0 

3.30 

31625 

741011 

5.00 

4.78 

-0.22 

60.2 

186.0 

33.0 

3.50 

87101 

741103 

3.10 

5.32 

0.22 

81.5 

181.0 

33.0 

3.50 

40240 

741203 

5.20 

6.39 

0.19 

72.5 

18  3-0 

33.0 

3.60 

67970 

741216 

5.10 

5.23 

0.  13 

71.8 

195.0 

33.0 

3.60 

SK3  P48E  IS  BEST  QUALITY  PBACTICABI* 
OJU-v  «o«JLLattJ»U  TO 


LASA  MAGNITUDE,  DMR=NI.ASA-f»  NETS,  SECTOR  L07 

24  EVENTS  APE  ANALYZED  FROM  210.1  TO  250.0  AZI 


EPX  NO. 

DATE 

MNEIS 

MLASA 

ML-MN 

DTST 

AZI 

DEPTH 

BFACTOR 

45075 

740123 

5.40 

5.39 

-0. 01 

95.6 

242.0 

449.0 

3.60 

79770 

740228 

5.30 

4.60 

-0.  70 

89.8 

241.0 

98.0 

3.60 

340 

7t  ,1323 

5.30 

5.6? 

0.  33 

91.5 

240.0 

33.0 

3.80 

81970 

740517 

5.30 

4.92 

-0.38 

98.8 

240.0 

494.0 

3.90 

3235 

34  06  0 2 

5.30 

4.88 

-0.45 

87.1 

241.0 

7.0 

3.70 

3957 

740603 

5.40 

5.  81 

0.41 

93.0 

242.0 

33.0 

3.90 

33160 

740620 

5.30 

4.  98 

-0.32 

98.9 

241.0 

540.0 

3.90 

33860 

740620 

5.40 

5.  1° 

-0.21 

36.3 

243.0 

31.0 

3.60 

46135 

740626 

5.40 

5. 09 

-0.  31 

97.3 

242.0 

551.0 

3.70 

69470 

740714 

5.. 30 

4.25 

-1.05 

99.0 

237.0 

183.0 

3.90 

77160 

740710 

5.30 

4.85 

-0.45 

e6.i 

241.0 

33.0 

3.60 

81715 

740723 

5.50 

5.35 

-0.  15 

87.4 

242.0 

46.0 

3.70 

83210 

740724 

5.40 

5.20 

-0.20 

101.  1 

235.0 

33.0 

4.10 

15030 

740810 

5.50 

5.37 

-0.13 

94.  4 

243.0 

602.0 

3.60 

4"M85 

740825 

5.30 

4.95 

-0.35 

96.5 

242.0 

542.0 

3.70 

52615 

74082° 

5.  30 

5.  10 

-0.20 

99.  1 

237.0 

58.0 

4.10 

03135 

740920 

5.30 

5.38 

0.08 

°4. 3 

239.0 

33.0 

3.80 

38830 

741014 

5.30 

3.87 

-1.43 

93.3 

242.0 

554.0 

3.50 

49140 

741018 

5.40 

5.28 

-0.12 

86.5 

241.0 

33.0 

3.70 

95510 

741108 

5.30 

5.21 

-0.09 

86.4 

242.0 

12.0 

3.60 

96520 

741108 

5.40 

8.56 

0.16 

86.4 

242.0 

3 3.0 

3.60 

23430 

741123 

5.40 

4.96 

-0.44 

94.5 

241.0 

211.0 

3.70 

°3875 

740729 

5.40 

5.76 

0.36 

91.4 

245.0 

586.0 

3.50 

21950 

740813 

5.40 

5.03 

-0.37 

90.5 

247.0 

56.0 

3.70 

* 1 . Jl 

mis  ftofeafc  bust  qualI  rr  practicabij 

• FRO*  COPY  rURtiiSKSC  ro  0DQ  - — 

t» 

t LASA  A AONTTHDE,  D*  8=  ALA  S A-ANFIS  , S EC^OR  L08 

« 17  EVENTS  APE  ANALYZED  FROM  250.1  TO  290.0  AZI 


EPX  NO. 

DATE 

ANEIS 

ALAS  A 

AL-AN 

DIST 

AZI 

DEPTH 

BFACTOR 

43296 

740121 

5.40 

5.52 

0.  12 

93.5 

266.0 

38.0 

4.20 

54240 

740201 

5.40 

5.05 

-0.35 

100.8 

270.0 

33.0 

4.00 

55400 

740201 

5.40 

4.92 

-0.48 

101.0 

271.0 

47.0 

4.00 

57355 

740203 

5.40 

4.91 

-0.49 

101.0 

270.0 

43.0 

4.00 

59  365 

740204 

5.40 

4.73 

-0.67 

100.7 

270.0 

55.0 

4.10 

63060 

740210 

5.40 

4.85 

-0.55 

101.2 

271.0 

3 8.0 

4.00 

9845 

740402 

5.40 

5.20 

-0.20 

100.8 

271.0 

47.0 

4.10 

58995 

740504 

5.50 

5.44 

-0.06 

94.1 

254.0 

602.0 

3.60 

12470 

740608 

5.40 

5.  50 

0.  10 

99.0 

265.0 

34  ^0 

4.20 

18140 

740612 

5.40 

5.25 

-0.15 

99.3 

265.0 

33.0 

4.20 

95160 

740725 

5.50 

5.19 

-0.  31 

101.7 

273.0 

33.0 

4.10 

22185 

740813 

5.50 

5.41 

-0.09 

102.7 

275.0 

100.0 

4.30 

25340 

740815 

5.50 

5.63 

0.  13 

99.7 

267.0 

59.0 

4.10 

1055 

740925 

5.40 

5.45 

0.05 

106.5 

273.0 

10.0 

4.60 

14895 

741119 

5.50 

4.60 

-0.90 

101.3 

277.0 

18.0 

4.00 

49650 

741207 

5.40 

4.96 

-0.44 

100.5 

272.0 

93.0 

4.  10 

91080 

741231 

5.50 

6.07 

0.  57 

49.6 

252. C 

5.0 

3.50 

r 

f 


XBX3  PAGE  IS  BEST  QUALITY  PRACTICABLI 
efts'!  lutuu-oiiJw  TO  DOC  


LASA  MAGNITUDE,  DHB=t!LASA-MNEIS , SECTOR  L09 

151  EVENTS  ARE  ANALYZED  FPOM  290.1  TO  325.0  AZI 


EPX  NO. 

DATE 

NNEIS 

MLASA 

ML-MN 

DIST 

AZI 

DEPTH 

BFACTOR 

37915 

740622 

5.20 

4.74 

-0.46 

73.8 

311.0 

40.0 

3.50 

33740 

740622 

5.10 

5.  11 

0.01 

87.9 

295.0 

103.0 

3.50 

46370 

740627 

*.00 

4.  36 

-0.62 

85.  1 

313.0 

39.0 

3.70 

59655 

740707 

5.20 

5.24 

0.04 

64.0 

311.0 

33.0 

3.70 

66645 

740713 

5.10 

5.11 

0.01 

93.0 

315.0 

97.0 

3.80 

“0050 

740728 

5.00 

5.14 

0.14 

63.8 

310.0 

60.0 

3.40 

95215 

740730 

5.00 

4.89 

-0.11 

64.0 

310.0 

42.0 

3.50 

98160 

740801 

5.20 

4.74 

-0.46 

29.6 

306.0 

33.0 

3.30 

99285 

740801 

5.00 

4.63 

-0.37 

87.6 

299.0 

321.0 

3.40 

46840 

740825 

5.00 

5.01 

0.01 

80.0 

306.0 

33.0 

3.40 

50685 

740828 

5.00 

4.74 

-0.  26 

77.3 

310.0 

51.0 

3.50 

56800 

740830 

5.30 

5.03 

-0.27 

81.3 

306.0 

24.0 

3.50 

56195 

740930 

5.10 

4.27 

-0.83 

81.4 

306.0 

46.0 

3.40 

61845 

740903 

5.30 

4.98 

-0.32 

79.9 

306.0 

20.0 

3.40 

63440 

740903 

5.00 

4.60 

-0.40 

79.9 

306.0 

15.0 

3.40 

65890 

740905 

5.10 

4. 97 

-0.13 

79.8 

306.0 

33.0 

3.40 

71630 

740908 

5.10 

5.27 

0.  17 

86.9 

299.0 

135.0 

3.50 

8*7  35 

740916 

5.30 

5.01 

-0.29 

67.6 

311.0 

54.0 

3.50 

86010 

740917 

5.00 

4.  50 

-0.50 

29.2 

306.0 

17.0 

3.30 

16505 

741004 

5.00 

4.45 

-0.55 

56.7 

313.0 

33.0 

3.50 

33795 

741012 

5.30 

5.  32 

0.02 

73.0 

311.0 

26.0 

3.60 

39835 

741015 

5.40 

4.98 

-0.42 

72.8 

311.0 

22.0 

3.60 

40240 

741015 

5.00 

4.  80 

-0.20 

67.7 

311.0 

79.0 

3.30 

64515 

741024 

5.10 

4.20 

-0.90 

81.3 

306.0 

45.0 

3.40 

79700 

741030 

5.  30 

5.  50 

0.20 

88.0 

313.0 

33.0 

3.80 

79885 

741030 

5.00 

4. 86 

-0.14 

85.3 

304.0 

422.0 

3.30 

84725 

741102 

5.10 

4.  60 

-0.50 

73.0 

313.0 

67.0 

3.30 

85810 

741102 

5.10 

4.54 

-0.56 

77.1 

310.0 

44.0 

3.60 

92950 

741107 

5.00 

5.42 

0.42 

69.5 

312.0 

55.0 

3.50 

97995 

741109 

5.40 

4.69 

-0.71 

73.4 

311.0 

27.0 

3.60 

670 

741111 

5.00 

5.05 

0.05 

89.6 

313.0 

44.0 

3.70 

12180 

741117 

5.00 

4.52 

-0.48 

72.5 

312.0 

42.0 

3.60 

30795 

741128 

5.30 

4.75 

-0.55 

36.4 

302.0 

32.0 

3.30 

43705 

741205 

5.00 

4.58 

-0.  42 

89.4 

295.0 

185.0 

3.30 

45390 

741205 

5.00 

4.91 

-0.09 

81.0 

305.0 

3 3.0 

3.50 

55730 

741210 

5.00 

4.98 

-0.02 

70.0 

311.0 

36.0 

3.60 

64200 

741213 

5.00 

4.58 

-0.42 

62.9 

311.0 

33.0 

3.60 

76845 

741225 

5.10 

5.51 

0.41 

49.4 

306.0 

37.0 

3.50 

89125 

741231 

5.10 

5.31 

0.21 

45.6 

305.0 

228.0 

2.90 

62885 

74  05  07 

5.00 

4.40 

-0.60 

96.6 

316.0 

36.0 

4. 10 

59905 

740505 

5.50 

4.85 

-0.6* 

98.7 

316.0 

26.0 

4.20 

25383 

740U12 

5.50 

*.uo 

-0.01 

98.5 

301.0 

33.0 

4.20 

38615 

740116 

5.10 

4.24 

-0.86 

97.7 

316.0 

58.0 

4.00 

56255 

740502 

5.00 

4.91 

-0.09 

78.2 

309.0 

33.0 

3.60 

16970 

741004 

5.  30 

4.56 

-0.74 

56.6 

312.0 

16.0 

3.50 

97410 

740801 

5.20 

4.49 

-0.71 

29.6 

306.0 

10.0 

3.30 

8810 

740402 

5.20 

5.26 

0.06 

73.4 

313.0 

41.0 

3.60 

97905 

740322 

5.10 

5. 13 

0.03 

36.2 

302.0 

33.0 

3.30 

6680 

740806 

5.00 

4.52 

-0.48 

30.3 

31  3.0 

136.0 

3.20 

24105 

741123 

5.40 

5.25 

-0.15 

96.5 

315.0 

33.0 

4.10 

41620 

740823 

5.00 

4.75 

-0.25 

97.3 

316.0 

33.0 

4.10 

12615 

740808 

5.40 

4.92 

-0.48 

96.1 

316.0 

32.0 

4.00 

26215 

740617 

5.10 

5.05 

-0.05 

6 2.  1 

311.0 

33.0 

3.70 

4«280 

740127 

5.30 

5.18 

-0.12 

87.  1 

298.0 

148.0 

3.40 

61005 

740505 

5.10 

4.99 

-0.11 

66.  1 

312.0 

143.0 

3.30 

30955 

740105 

5.40 

5.40 

0.00 

41.  3 

302.0 

41.0 

3.20 

77205 

740415 

5.00 

4.73 

-0.27 

39.7 

302.0 

44.0 

3.10 

38755 

741014 

5.30 

4.45 

-0.85 

72.9 

311.0 

15.0 

3.60 

54290 

741020 

5.30 

5.40 

0.10 

72.3 

313.0 

24.0 

3.60 

64260 

741024 

5.40 

5.42 

0.02 

79.8 

11S:8 

63.0 

3.40 

71140 

741027 

5.  10 

5.  34 

0.24 

69. 1 

223.0 

3.20 

A— 10 


1 


* .*  V'V\  * 

*HIS  PAGE  Ifi'-MST'  QUALITY  PRACtLCABLH 
IRON  COPY  PURltU-SiM)  JO  DflQ  ^ 


LASA  NAGNTTtJDE,  OH D=4LA  SA -l» NETS  , SECTOR  L09 

151  EVENTS  AFF  ANALYZED  FPOH  290.1  TO  325.0  AZI 


EPY  NO. 

DATE 

HNEIS 

HLASA 

Hl-NN 

DIST 

AZI 

DEPTH 

BFACTOR 

5060 

741114 

5.50 

5.86 

0.36 

30.9 

310.0 

37.0 

3.40 

12460 

741117 

5.30 

5.17 

-0.13 

64.6 

315.0 

31.0 

3.50 

16660 

741120 

5.00 

5.23 

0.23 

37.3 

302.0 

57.0 

3. 30 

29675 

741120 

5.20 

0.61 

0.41 

43.6 

303.0 

63.0 

3.20 

34345 

741130 

5.20 

5.42 

0.22 

49.5 

309.0 

17.0 

3.40 

21075 

741007 

5.20 

l-M 

0.  32 

69.2 

315.0 

325.0 

3.20 

12275 

85805 

741002 

5.20 

0.19 

58. 1 

1U:S 

51.0 

3.50 

740916 

5.50 

5.95 

0.45 

60.4 

48.0 

3.50 

61916 

740903 

5.20 

5.70 

0.50 

60.3 

312.0 

306.0 

22.0 

40.0 

3.50 

46245 

740025 

5.30 

5.29 

-0.01 

PO.O 

3.40 

32960 

740810 

5.00 

5.42 

0.42 

49.0 

305.0 

33.0 

3.  40 

55 

740801 

5.30 

c . 73 

0.43 

60.  3 

312.0 

ttl.O 

3.60 

57175 

741230 

5.10 

0.37 

0.27 

20.8 

317.0 

62.0 

3.10 

34000 

7u 1 0 1 2 

5.50 

5.46 

-0.06 

73.0 

311.0 

24.0 

3.60 

40650 

741203 

6.00 

0.32 

0.32 

62.6 

312.0 

36.0 

3.60 

64740 

740904 

5.30 

4.93 

-0.37 

74.  1 

312.0 

55.0 

.3.40 

46735 

740025 

5.20 

0.06 

-0.  16 

PO.O 

306.0 

81.0 

3.  30 

45705 

740024 

5.70 

4.64 

-1.06 

99.6 

301.0 

37.0 

.3.10 

37060 

7404  2 1 

5.30 

5.56 

0.26 

67.8 

314.0 

28.0 

3.70 

3^710 

740610 

5.00 

5.51 

0.  51 

54.  1 

31  3.0 

46.0 

3.50 

36580 

740R19 

5.30 

5.38 

0.08 

00.5 

309.0 

23.0 

3.50 

23515 

740014 

5.30 

0.  40 

0.  19 

72.  9 

312.0 

41.0 

3.60 

76°  6 

74  08  06 

5.20 

6.81 

-0.39 

65.0 

312.0 

166.0 

3.20 

23985 

7404  1 1 

5.  30 

5.  12 

-0.  18 

71.  1 

312.0 

75.0 

3.30 

20910 

740409 

6.20 

6.36 

0.14 

06.2 

295.0 

U3.0 

3.50 

.20785 

740409 

6.80 

5.77 

0.27 

67.0 

312.0 

139.0 

3.40 

12020 

7^0404 

5.30 

5.28 

-0.02 

76.  c 

311  .0 

97.0 

3.  30 

45315 

7 U 0 3 1 9 

5.00 

4.96 

-0.06 

78.  4 

30Q.0 

69.0 

3.30 

4630 

740806 

0.18 

5.36 

0.24 

42.4 

30  3.0 

66.0 

3.30 

92390 

740315 

0.00 

5.  4? 

0.42 

r 9. 2 

311.0 

35.0 

3.50 

86840 

740300 

5.00 

5.70 

0.26 

06.5 

313.0 

42.0 

3.70 

74915 

7408 ^ 3 

0.30 

5.07 

-0.21 

72.0 

313.0 

64.  C 

3.30 

8 14  5 5 

740778 

5.00 

0.31 

0.31 

63.0 

311.0 

54.0 

3.50 

80880 

740720 

6.  10 

5.37 

0.  27 

64.0 

310.0 

52.0 

3.50 

2U7c0 

741008 

5.00 

4.5P 

-0.42 

78.2 

311.0 

56.0 

3.  30 

481^0 

740126 

*.30 

5.50 

0.20 

01.2 

102.0 

c 4 . 0 

3.10 

41590 

740118 

5. On 

5.12 

0.12 

39.  1 

302.0 

5 9.0 

3.  30 

62015 

740203 

5.40 

5.  36 

-0.02 

C1.P 

312.0 

33.0 

7.40 

M 26  0 

7402  07 

0.00 

5.56 

0.56 

59.9 
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APPENDIX  B 

Signal  Variance  Among  Sources  and  Receivers  as  a 
Consequence  of  the  Seismic  Reciprocity  Theorem 
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From  Knopoff  and  Gangi  (1959),  the  seismic  reciprocity  theorem  for  linear 
elasticity  is 

l ui(P,Q)  ft(P)  = Uj(Q,P)  f j (Q) 

where  u^(P,Q)  is  the  i component  of  displacement  at  receiver  point  P due  to 
source  at  Q of  strength  f ^ (Q)  in  the  direction  j.  The  theorem  is  valid  in 
inhomogeneous,  non-isotropic  media  with  arbitrary  boundaries. 

In  this  form  the  theorem  is  not  suited  for  applications  to  problems  of 
magnitude-yield  because  an  explosion  source  is  not  a directed  force,  but  rather 
a dilatation  expanding  against  hydrostatic  pressure.  The  relevant  generali- 
zation is  (Knopoff,  1978  personal  communication) 

• • • 

V»u(P,Q)  t(P)  = V«u(Q,P)  t(Q) 

Verbally,  this  equation  can  be  expressed  as:  a unit  dilatation  at  Q work- 
ing against  hydrostatic  pressure  x(Q)  produces  a divergence  at  P that  when 
multiplied  by  t(P)  is  equal  to  the  divergence  produced  at  Q by  a unit  dila- 
tation at  P multiplied  by  t(Q). 

4 

Again,  this  theorem  is  valid  in  inhomogenous , non-isotropic  media  with 
I arbitrary  boundaries. 

Explosions  of  equal  yield  at  equal  depth  should  satisfy  the  condition  of 
equal  dilatations  at  equal  pressures.  Thus  in  applying  this  to  the  present 
study,  if  the  divergence  (of  a signal  from  Semipalatinsk)  could  be  measured 
at  a shot  d^pth  at  subarrays  of  LASA;  the  variance  of  the  divergence  measure- 
ments would  be  equal  to  the  variance  observed  in  the  divergence  measured  at 
shot  depth  near  the  shot  at  Semipalatinsk  from  a series  of  explosions  set  off 
at  the  subarrays  of  LASA. 

The  only  problem  here  is  that  in  practice  we  measure  displacements  at  the 
surface,  not  divergence  at  shot  depth.  However,  to  the  extent  that  the  incom- 
ing wave  at  both  LASA  and  Semipalatinsk  can  be  modelled  as  a plane  wave  the 
divergence  may  be  simply  calculated  as  the  vertical  derivative  of  the  incoming 
t wave  amplitude  with  a delayed  signal  of  the  same  polarity  (the  reflection 

from  the  surface).  If  focussing  and  de-focussing  are  the  main  causes  of  signal 
* amplitude  variation,  then  the  maximum  divergence  at  depth  will  be  a constant 
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factor  times  the  maximum  displacement  at  the  surface.  Thus,  the  variance  of 
the  divergence  in  logarithmic  units  will  be  the  same  as  the  variance  of  the 
displacement.  Thus,  finally,  the  variance  of  the  displacement  at  Semipalatinsk 
from  shots  spaced  around  LASA  will  be  equal  to  the  variance  observed  at  the 
shotpoints  from  a shot  at  Semipalatinsk. 

In  this  report  variance  measurements  have  been  made  at  LASA  on  the  signals 
from  earthquakes,  not  explosions.  Thus,  while  the  theorem  is  not  totally 
valid,  we  can  see  that  our  conclusions  about  variance  are  still  valid.  Consider 
that  the  P-wave  signal,  as  it  emerges  from  the  earthquake,  could  have  come  from 
an  explosion  of  suitably  determined  yield.  In  logarithmic  units  this  factor 
would  cancel  out  of  variance  estimates  so  that  we  could  predict  the  variance 
to  be  observed  at  the  earthquake  foci  from  shots  spaced  around  LASA.  Since 
the  crust  above  the  earthquake,  between  the  earthquake  and  a surface  sensor,  is 
a complex  part  of  the  path,  the  variance  observed  at  the  surface  would  then  be 
greater  than,  or  equal  to,  the  variance  observed  at  LASA  from  the  earthquakes. 
This  conclusion  applies  only  so  long  as  the  radiation  patterns  of  the  earth- 
quakes do  not  vary  significantly  across  LASA. 

Note  also  that  the  variance  observed  at  RKON  due  to  shots  spaced  around 
the  Nevada  Test  Site  (NTS)  is  equal  to  the  variance  that  would  be  observed  at 
NTS  from  a shot  at  RKON.  Or,  another  example,  the  variance  observed  at  NTS 
from  a shot  at  Semipalatinsk  is  equal  to  the  variance  which  would  be  observed 
at  Semipalatinsk  from  shots  spaced  around  NTS.  Clearly,  there  are  many  appli- 
cations of  this  theorem. 


